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I.  FORWARD 


NT  IS  CRA&I 
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A  major  objective  of  our  research  under  the  ARO  contract  is  to  investigate  novel  ap- T 
proaches  which  will  lead  to  practical,  high  frequency  (>  100  GHz)  sources.  For  cost  effec¬ 
tiveness,  we  have  carried  out  a  significant  portion  of  the  work  at  lower  frequencies  where  we 
can  economically  assemble  the  basic  experiments  necessary  to  carry  out  proof-of-principle 
physics  studies.  One  must  demonstrate  ingenuity  and  cleverness  to  make  dramatic  advances 
with  our  moderate  level  of  funding.  However,  as  discussed  in  detail  in  Secs.  Ill  and  IV.  all 
of  the  devices  under  study  have  the  potential  for  high  frequency  operation. 

During  the  last  year  of  the  contract,  four  important  papers  describing  our  ARO-funded 
contributions  to  millimeter- wave  and  sub-millimeter  wave  source  development  were  pub¬ 
lished  in  Special  Issues  of  Scholarly  Journals.  The  paper,  “Marginal  Stability  Design  Crite¬ 
rion  for  Gyro-TWTs  and  Comparison  of  Fundamental  with  Second  Harmonic  Operation,” 
by  A.T.  Lin,  K.R.  Chu,  C.C.  Lin,  C.S.  Kou,  D.B,  McDermott  and  N.C.  Luhmann,  Jr., 
{Int.  J.  Electronics,  72,  873  (1992))  introduced  a  method  for  building  stable  gyro-TWTs 
with  the  best  performance  possible.  The  papers,  “Iligh-Power  Harmonic  Gyro-TWT’s  - 
Part  I:  Linear  Theory  and  Oscillation  Study,”  by  C.S  Kou,  Q.S.  Wang,  D.B.  McDermott. 
A.T.  Lin,  K.R.  Chu,  and  N.C.  Luhmann,  Jr.,  (IEEE.  Trans  on  Plasma  Science,  20,  155, 
(1992))  and  “High-Power  Harmonic  Gyro-TWTs  -Part  II:  Nonlinear  Theory  and  Design.” 
by  Q.S.  W'ang.  C.S.  Kou.  D.B.  McDermott.  .A.T.  Lin.  K.R.  Chu.  and  N.C.  Luhmann.  Jr. 
(IEEE  Trans,  on  Plasma  Science.  20,  163  (1992)),  explained  why  harmonic  gyro-TWTs 
can  be  expected  to  deliver  extremely  high  power  and  yet  remain  stable.  Finally,  the  paper, 
“Bragg  Reflectors,”  by  C.K.  Chong,  D.B.  McDermott,  M.M.  Razeghi,  N.C.  Luhmann,  Jr., 
J.  Pretterebner,  D.  Wagner,  M.  Thumm,  B.  Kulke  and  M.  Caplan  (IEEE  Trans,  on  Plasma 
Science  2Q,  393  (1992)),  describes  the  performance  and  capability  of  a  frequency  selective 
reflector  that  is  sure  to  play  an  important  role  in  future  high  power  overmoded  oscillators. 

The  organization  of  this  final  report  is  as  follows.  Section  11  contains  a  list  of  the  report's 
illustrations,  tables  and  appendixes.  Section  HI  slates  the  objectives  of  the  research.  .A 
summary  of  the  research  findings  is  presented  in  Sec.  IV  with  references  given  in  Sec.  I\’.F'. 
Our  publications  and  technical  reports  produced  during  the  last  year  are  listed  in  Sec.  V. 
Section  \T  lists  all  scientific  personnel  at  UCLA  participating  in  the  research.  Section  \’I1 
contains  the  report  on  inventions.  Copies  of  four  publications  funded  by  the  contract  are 
included  in  the  Appendix. 
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III.  STATEMENT  OF  PROBLEM 


The  overall  objective  of  our  research  lor  ARO  is  to  develop  devices  that  have  the  poten¬ 
tial  for  being  compact  generators  of  submillimeter-wave  radiation.  Our  efforts  have  focussed 
on  gyrotron  devices  because  they  have  proven  to  be  efficient  generators.  High  frequency 
gyrotrons  can  be  made  more  compact  by  reducing  the  cathode  voltage  and  by  operating  at 
harmonics  of  the  cyclotron  frequency  which  reduces  the  required  magnetic  field  by  a  factor 
of  the  harmonic  number.  In  the  first  eight  years  of  our  .4R0  funding,  we  developed  high 
harmonic  gyro  devices  driven  by  rf- accelerated  electron  beams,  with  an  eleventh-harmonic 
gyrotron.  a  fifth-harmonic  gyro-klystron  amplifier  and  an  eighth-harmonic  gyro-TWT  am¬ 
plifier  being  three  notable  accomplishments.  It  is  of  interest  to  note  that  industry  (Varian, 
Litton  and  Northrop)  is  currently  developing  several  types  of  high-harmonic  gyro-amplifiers 
at  95  GHz,  based  in  large  part  on  our  work.  Additional  research  efforts  include  a  successful 
250  GHz  CARM  oscillator  experimental  coUaboration  with  Lawrence  Livermore  National 
Laboratory. 

During  the  most  recent  contract  period,  we  have  continued  our  development  of  gyro- 
devices  that  have  a  potential  for  high  frequency,  high  harmonic  operation.  (1)  The  high 
harmonic  gyro-frequency  multiplier  is  an  efficient  generator  of  high  frequency  rf.  where  the 
harmonic  number  can  be  as  high  as  ten.  or  perhaps  even  higher.  The  primary  problem  is 
to  inject  an  electron  beam  into  the  accelerator  with  a  beam  radius  much  smaller  than  the 
transverse  wavelength  of  the  mode  being  excited.  Another  major  problem  is  to  taper  the 
magnet  field  so  that  the  electrons  continuously  lose  energy  during  their  transit  through  the 
cavity.  (2)  Using  the  marginal  stability  design  criterion,  high-power  harmonic  gyro-TWTs 
can  be  designed  to  operate  in  high  order  modes  in  large  size  waveguide.  The  primary 
problem  is  to  determine  the  start-oscillation  lengths  for  all  modes  in  the  tube.  The  other 
major  problem  is  to  find  whether  there  is  sufficient  gain  for  the  gyro-TWT's  operating 
mode.  (3)  In  our  high  power,  second-harmonic  TE21  gyro-TWT.  the  primary  problem  will 
be  to  create  an  electron  beam  with  a  small  spread  in  its  axial  velocity.  Another  major 
problem  will  be  to  suppress  harmonic  gyro-BWO  modes  whose  instability  can  be  excited 
by  poorly  terminated  .severs.  (4)  The  high  harmonic  C.AR.M  will  also  be  dependent  on  a 
high  quality  electron  beam.  In  addition,  the  circuit  must  be  sliced  carefully  to  suppress  the 
competing  modes  without  interfering  with  the  design  mode.  (5)  In  the  initial  gyro-klystron 
experiment,  the  primary  problem  will  be  to  suppress  gyrotron  oscillation  in  the  cavities. 
The  next  major  problem  will  be  to  keep  the  beam  sufficiently  cold  so  that  high  efficiency 
can  be  realized.  In  the  succeeding  harmonic  gyro-klystron  experiments,  it  will  be  even  more 
challenging  to  suppress  the  excitation  of  the  stronger,  lower  harmonic  gyrotron  modes  in 
the  cavities  and  in  the  drift-tubes. 
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IV.  SUMMARY  OF  RESEARCH  FINDINGS,  6/25/91  -  6/24/92 


A.  High  Harmonic  Gyro- Frequency  Multiplier 
Summary 


A  major  objective  of  the  program  is  to  upgrade  our  previous  third-harmonic  gyro- 
frequency  multiplier  experiment  to  much  higher  output  power  and  efficiency  and  at 
a  much  larger  harmonic.  A  simulation  code  has  been  w'ritten  to  design  the  output 
cavity  and  investigate  efficiency  enhancement  through  magnetic  tapering.  Efficiencies 
as  high  as  70%  are  realizable.  A  new  buncher/accelerator  cavity  has  been  built  to 
accelerate  much  higher  levels  of  electron  beam  current.  A  new  high-perveance  electron 
beam  injector  has  been  built  and  successfully  tested  for  injecting  the  required  higher 
values  of  beam  current. 


We  have  been  involved  in  the  design  of  a  significant  extension  of  our  successful  high 
harmonic  gyro-frequency  multiplier  (i.e.,  prebunched  gyrotron).  In  addition  to  potentially 
generating  rf  in  the  sub-mm  wave  region,  the  gyro-multiplier  promises  to  be  extremely 
efficient  because  the  electrons  from  the  rf-accelerator/buncher  cavity  enter  the  output  cavity 
already  bunched  and,  due  to  the  extended  gyro  interaction,  they  remain  in  near  resonance 
throughout  the  cavity.  This  results  in  a  first-order  interaction,  rather  than  the  weaker 
second-order  interaction  associated  with  the  conventional  gyrotron.  The  results  of  our 
first  prebunched  gyrotron  experiment  [1]  was  that  our  third  harmonic  multiplier  generated 
28  GHz  rf  with  a  conversion  power  of  6.7  kW  with  an  efficiency  of  13%.  In  this  initial 
experiment,  the  beam’s  radial  guiding  center  spread  was  fairly  large  so  that  all  the  electrons 
did  not  enter  the  multiplier  cavity  at  the  same  decelerating  phase  as  desired. 

The  high  frequency  potential  of  this  device  is  primarily  limited  only  by  the  availability  of 
high  power,  high  frequency  sources  to  drive  the  accelerator /buncher  cavity  [2].  For  example, 
a  tenth-harmonic  gyro-multiplier  utilizing  a  140  GHz  rf-accelerator  would  generate  high 
power  electromagnetic  waves  with  a  wavelength  of  0.21  mm  in  a  magnetic  field  of  50  kG. 
Gyrotrons  currently  exist  which  can  generate  cw  power  levels  of  0.5  MW  at  140  GHz  and 
higher  frequency  gyrotrons  are  under  development.  Therefore,  this  frequency  multiplication 
approach  offers  exciting  possibilities  for  submillimcler  wave  power  generation. 

The  gocils  of  the  new  experiment  are  1 )  to  use  a  much  more  tightly  focused  electron 
beam  so  that  the  entire  beam  enters  the  multiplier  cavity  closer  to  the  same  optimum 
phase,  2)  to  use  a  much  higher  value  of  beam  current  for  higher  power.  3)  to  operate  at 
much  higher  harmonics,  and  4)  to  design  the  cavity  and/or  use  magnetic  tapering  for  higher 
efficiency.  We  have  decided  to  perform  the  next  prebunched  gyrotron  experiment  at  the 
tenth  harmonic,  which  will  be  at  a  frequency  of  29  GHz  since  a  2.9  GHz,  1  MW  klystron 
amplifier  will  be  used.  However,  as  mentioned  above,  this  concept  scales  to  much  higher 
frequencies.  We  have  built  a  new  gyroresonanl  rf  accelerator  cavity.  Whereas  the  cavity 
used  in  our  previous  experiments  was  designed  to  accelerate  electrons  to  high  energy  (250- 
500  keV),  the  new  cavity  is  designed  to  accelerate  much  higher  current  beams  (2-5  A)  to  the 
more  moderate  energy  of  50-  100  keV’.  The  length  of  the  new  accelerator  is  roughly  one- half 
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of  the  previous  cavity.  Also,  lossy  Eccosorb  epoxy  has  been  added  to  lower  its  loaded  Q  to 
150.  The  cavity  is  currently  undergoing  cold-tests  to  measure  its  Q  factor  while  we  increase 
the  coupling  by  enlarging  the  size  of  the  coupling  slots.  Our  simulation  code  has  been  used 
to  study  the  acceleration  characteristics  of  the  new  cavity.  Figure  1  shows  the  performance 
that  can  be  expected  when  using  75%  of  our  available  input  power  of  1  MW  .  As  desired,  the 
final  electron  energy  and  velocity  ratio  o  (=  ri/i’||)  are  relatively  constant  as  the  magnetic 
field  is  varied  around  the  operating  value  of  Oco/^W  =  1-2. 


Figure  1;  Predicted  dependence  of  Lorentz  factor  and  velocity  ratio  o  =  t’i/nn  on  ratio  of 
rest-mass  cyclotron  frequency  to  accelerator  rf  frequency  in  gyroresonant  TEm  rf  acceler¬ 
ator.  The  field  value  of  A  =  0.03  would  be  achieved  for  Q=I35  with  a  power  of  P  =  750 
kW. 

We  have  devoted  considerable  effort  to  construct  a  high  perveance  electron  gun  and 
optimize  the  magnetic  pole  piece.  The  gun  is  identical  to  the  electron  gun  from  the  SL.^C 
5045  Klystron,  but  scaled  down  in  size  by  a  factor  of  3.5  for  lower  voltage  operation.  In 
order  to  achieve  the  same  magnetic  field  shape  as  in  the  5045  Klystron  (but  scaled  by  the 
same  factor  of  3.5),  our  magnetic  pole  piece  has  a  DC  offset  coil  and  a  Helmholtz  gradient 
coil  in  addition  to  the  iron  pole  piece.  As  shown  in  Fig.  2,  the  resultant  magnetic  field  shape 
is  a  very  close  fit  to  the  desired  profile,  especially  at  the  emitting  surface,  which  is  the  most 
important  region.  With  this  Pierce  electron  gun,  we  were  able  to  successfully  transmit  a 
5A,  30  kV  beam  through  a  3  mm  diameter  aperture.  In  addition,  the  solenoid  has  been 
reconfigured  to  contain  a  separate  solenoid  for  the  electron  gun  and  a  long  uniform  region 
for  the  gyroresonant  rf  accelerator  and  gyro-multiplier  cavities.  The  magnetic  field  profile 
of  the  new  solenoid  is  shown  in  Fig.  3. 

A  particle  tracing  simulation  code  was  written  to  design  an  optimized  frequency  multi¬ 
plier  cavity  for  the  new  experiment.  To  first  order,  where  the  finite  beam  size  is  neglected, 
all  electrons  follow  the  same  trajectories  in  the  multiplier  cavity,  and  the  structure  of  the 
cavity  mode  is  unaffected  by  the  beam.  After  the  amplitude  of  the  mode  is  specified,  the 
code  follows  the  electrons  through  the  cavity  by  solving  the  equations  of  motion.  In  this 
way.  a  plot  of  efficiency  versus  wave  amplitude  can  be  constructed,  which  can  then  be  made 
into  a  plot  of  efficiency  versus  beam  current.  This  is  a  similar  procedure  to  that  used  for 
gyrotron  optimization  [3] 

Figure  4  shows  one  such  plot  for  a  cavity  corresponding  roughly  to  our  previous  exper- 
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Axial  Position  (cm) 


Figure  2:  Dependence  of  magnetic  field  in  electron  gun  for  scaled  SLAC  5045  gun  from 
measurement  (solid  curve)  and  scaling  of  the  actual  SLAC  5045  profile  (dashed  curve).  The 
center  of  the  emitting  surface  is  at  z  =  8.8  cm. 


Axial  Position  (cm) 


Figure  3;  Profile  ol  magnetic  field  in  .solenoid  for  high-harmonic  gyro-frequency  multiplier 
experiment. 


Figure  4:  Dependence  of  conversion  efficiency  on  phase  of  TE312  wave  at  input  of  third 
harmonic  multiplier  for  ten  values  of  normalized  electric  field  amplitude,  £  =  Eoer^jmc^ 
(200  kV,  a  =  2.0,  L/a  =  2.0,  E  =  0.2,  0.4 . 2.2). 


iment.  If  the  electrons  enter  the  cavity  at  the  optimum  small-signal  wave  phase  of  (j>  =  rc 
(maximum  deceleration),  then  the  peak  efficiency  is  30%.  The  electrons  fall  out  of  resonance 
as  they  lose  energy.  An  efficiency  of  nearly  70%  is  obtained  if  the  electrons  enter  the  cavity 
at  the  phase  at  which  the  electric  field  is  a  null.  Although  this  situation  probably  is  not 
realizable  since  the  wave  will  grow  for  the  phase  which  gives  maximum  small  signal  gain 
(4)  =  tt),  a  similar  efficiency  will  result  if  the  magnetic  field  is  tapered.  Magnetic  tapering 
will  be  investigated.  The  purpose  of  making  these  plots  is  to  determine  the  cavity  Q  at 
which  the  optimum  efficiency  will  be  attained  for  the  available  beam  current.  The  cavity’s 
length  will  be  varied  to  find  the  value  for  optimum  efficiency.  We  will  also  investigate  the 
effectiveness  of  magnetic  field  tapering  to  prolong  resonance  and  thereby  enhance  the  effi¬ 
ciency.  Since  our  solenoid  is  comprised  of  many  discrete  magnets,  magnetic  tapering  is  an 
option.  Our  code  can  also  model  a  realistic  beam  with  finite  radial  thickness.  The  beam’s 
size  determines  the  maximum  harmonic  for  efficient  operation. 

Our  well-collimated  electron  beam  will  allow  efficient  operation  at  the  tenth  harmonic. 
Using  the  particle-tracing  simulation  code,  we  have  begun  to  design  a  29  GHz,  TEjo,i,2 
mode  tenth-harmonic  gyro-multiplier  cavity.  .Although  this  proof-of-principle  experiment 
will  generate  power  at  29  GHz  by  using  a  2.9  GHz  rf-accelerator/buncher.  it  must  be 
noted  that  power  at  much  higher  frequency  can  be  generated  by  using  a  higher  frequency 
accelerator/buncher  cavity.  For  example,  by  using  a  0.5  MW  cw,  140  GHz  gyrotron  to  drive 
the  accelerator/buncher  cavity,  a  scaled  ten-harmonic  multiplier  will  generate  cw  power  at 
1.4  THz. 
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B.  Overmoded  Gyro-TWT  Stability 


Summary 


A  gyro-TWT  amplifier  can  be  stabilized  by  breaking  it  into  several  short  stable  in¬ 
teraction  sections  separated  by  attenuating  severs  for  isolation.  The  marginal  stabil¬ 
ity  design  procedure  can  be  used  to  design  any  arbitrary  s'^-harmonic  TE,;,„-mode 
gyro-TWT  by  first  determining  the  start-oscillation  length  of  the  relevant  harmonic 
gyrotron  and  gyro-BWO  competing  modes.  The  length  of  the  interaction  sections 
must  be  less  than  the  shortest  start-oscillation  length.  The  number  of  interaction 
sections  is  then  chosen  to  obtain  the  desired  gain. 

High  power,  submillimeter-wave  devices  are  more  readily  built  if  the  interaction  waveg¬ 
uide’s  transverse  dimensions  are  as  large  as  possible.  To  increase  a  gyro-TWT’s  lateral 
dimensions,  it  is  necessary  to  increase  the  operating  mode  number  of  the  waveguide.  For 
example,  the  inner  diameter  of  a  110  GHz  gyro-TWT  designed  for  the  TE52  mode  would 
be  1  cm,  whereas  it  would  be  only  2  mm  if  the  lowest  order  TEn  mode  were  employed. 
A  conventional  device’s  stability,  which  is  problematic  even  in  a  fundamental  mode  de¬ 
vice,  would  be  highly  unlikely  for  such  an  overmoded  case.  However,  the  marginal  stability 
design  criterion  [4]  that  we  developed  considerably  relaxes  the  design  constraints.  Using 
this  procedure,  an  overmoded  gyro-TWT  would  be  comprised  of  many  short  stable  isolated 
sections. 

We  are  investigating  whether  the  marginal  stability  design  criterion  can  be  used  to  build 
a  high  performance,  fundamental  frequency  TEoi  gyro-TWT,  which  weis  the  configuration 
of  the  first  gyro-TWTs  to  be  tested.  [5]  If  they  can  be  kept  stable,  the  TEon  modes  are 
probably  the  best  modes  to  operate  a  gyro-TWT  because  of  their  low  ohmic  dissipation  and 
their  favorable  interaction  strength.  We  will  use  the  new  technique,  the  marginal  stability 
design  procedure,  to  design  a  stable  TEqi  gyro-TWT  amplifier.  For  the  TEoi  device,  the 
dominant  competing  mode  that  must  be  suppressed  is  the  TE21  gyro-BWO.  Therefore,  the 
interaction  sections  of  the  gyro-  TWT  must  be  chosen  shorter  than  the  start-oscillation 
length  for  the  TE21  gyro-BWO.  One  must  then  determine  how  much  gain  for  the  TEqi 
mode  will  result  for  that  length.  It  will  almost  certainly  be  less  than  the  gain  required  for 
most  applications  (>  30  dB).  Higher  gain  can  be  achieved  by  using  several  stable  sections 
separated  by  attenuating  severs. 


C.  High  Power,  Second  Harmonic  Gyro-TWT 


Summary 


An  extremely  high  power  three-section  second-harmonic  35  GHz  gyro-TWT  has  been 
designed  with  a  self-consistent  nonlinear  simulation  code  [6].  ,4n  output  power  of 
530  kW  with  21  %  efficiency,  54  dB  gain  and  6  %  bandwidth  is  predicted.  To  substan¬ 
tiate  the  theoretical  predictions  that  a  harmonic  gyro-TWT  can  generate  significantly 
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higher  output  power  with  better  stability  than  fundamental  gyro-TWTs  due  to  the 
higher  beam  currents  allowed  for  the  weaker  harmonic  interaction,  a  proof-of-principle 
Ku-band  experiment  is  being  constructed  at  UCLA  with  comparable  performance  pre¬ 
dicted.  A  single  anode  100  kV.  20  A  MIG  is  being  built  to  generate  the  a  =  t’i/e||  =  1 
electron  beam  with  At;||/i’||  =  8%.  The  input  and  rf  couplers  have  been  designed  and 
are  currently  being  built. 

Operation  at  the  second  harmonic  allows  the  requisite  magnetic  field  to  be  reduced  by 
a  factor  of  two,  an  important  concern  for  high  frequency  gyro-TWTs.  Combining  this  with 
the  case  used  in  the  previous  section  could  result  in  a  220  GHz,  TE52  gyro-TWT  amplifier 
with  an  inner  diameter  of  0.5  cm  in  a  magnetic  field  of  ~  45  kG.  Also,  it  must  be  pointed 
out  that  several  years  ago  we  successfully  tested  an  eighth-harmonic  gyro-  TWT  [7]  under 
ARO  funding. 

Harmonic  gyro-TWT’s  not  only  require  a  substantially  weaker  magnetic  field  and  offer 
nearly  the  same  high  efficiency  as  conventional  fundamental  gyro-i  .VT’s,  but  also  have  the 
potential  to  stably  generate  significantly  higher  levels  of  power  [4,  6,  8].  Basically,  since  the 
harmonic  interaction  is  generally  weaker  than  the  fundamental  interaction,  the  threshold 
beam  current  level  for  oscillation  is  raised  significantly.  This  is  seen  in  Fig.  5  where  the 
threshold  beam  current  for  the  absolute  instability  is  given  for  both  a  fundamental  and  a 
second-harmonic  interaction  (Bg  is  the  magnetic  field  at  which  the  particle  wave  resonance 
fine  grazes  the  operating  waveguide  mode  curve).  It  is  seen  that  for  o  =  v^/v\\  =  1,  the 
allowable  electron  beam  current  in  a  second  harmonic  gyro-TWT  is  higher  by  more  than 
fivefold  from  that  in  a  fundamental  gyro-TWT.  The  higher  working  beam  current  not  only 
provides  a  larger  amount  of  beam  power  available  for  millimeter-wave  power  generation, 
but  also  makes  up  for  the  loss  in  interaction  strength.  Thus,  a  high  performance  amplifier 
with  superior  power  capability  is  possible. 


Figure  5:  Dependence  of  threshold  beam  current  for  absolute  instability  on  velocity  ratio 
a  =  ux/u||  in  (a)  fundamental  TE\i  and  (b)  second-harmonic  T E21  gyro-TWT  amplifiers 
for  two  values  of  magnetic  field  (100  kV.  axis-encircling  beam). 
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Th*>  colt!  w.ivcpuitif  dispersion  r.'lation  of  tin-  7  pyro  TW'I'  aiiijdilicr  {llit*  siijiersrriin 
refers  tf»  the  cvclittutn  lianuonir  tiiinii»fr}  slum  it  in  (>  rt*v(*als  t  iie  key  feature  that 
inierat  t:on  at  tht*  fundaiuenia!  cyclotron  fn’»|iu*ncy  will  not  occur.  Tlt<*  u|itiiiiiini  olectroii 
beam  for  a  Rvro-lW  F  is  an  axis-enciicliiif'  electron  beam,  becau.se  it  would  produce 


the  strtuige.M  nueraction  anti  only  /  interactions  \rn  >  2}  coultl  then  occur.  Howt-ver, 
since  there  has  not  yet  been  a  successful  deinonstration  of  such  a  beam  with  the  requisite 
low  avial  velocity  spread,  the  liesipiis  f;iveii  in  this  paper  utilize  .\1I(1  Iteams  which  have 
long  been  employed  in  the  gyrotron  community  with  great  success. 


W 


Figure  G:  Uncoupled  dispersion  relation  of  operating  mode  (intersection  of  unbroken  curves) 


and  likely  oscillating  modes  (intersections  of  broken  curves  w’ith  negative  A.||)  for  TEji  6>'ro- 
TWT  amplifier  ( 100  kV.  o  =  1,  Eq  =  0,98Bg). 


Based  on  a  marginal  stability  procedure  (-I],  a  high  power  second  harmonic  TE^i  gyro- 

TWT  using  a  100  kV,  25  A  MIG  beam  with  a  =  1  and  At'||/e||  =  5  %  has  been  designed  (Gj, 

where  the  interaction  tube  is  comprised  of  three  sections  separated  by  absorbing  severs  (9). 

The  length  of  each  interaction  section  is  shorter  than  the  length  .at  which  the  device  will 

spontaneously  oscillate  and  is  determined  through  a  combination  of  analytical  theory  and 

multi-mode  PIC  code  simulation  [G].  Figure  7  shows  llie  threshold  lengths  for  the  most 

#2) 

threatening  oscillation  modes  in  our  T Ej,'  amplifier  design.  Here,  the  length  of  the  interac¬ 
tion  sections  is  limited  by  the  start  oscillation  Icngtli  for  the  3rrf  harmonic  TEjj  gyro-BWO 
mode.  The  number  of  .sections  is  chosen  to  yield  the  desired  gain.  For  the  530  kW  TEi, 
gyro-TWT  described  in  Table  1.  the  .spatial  power  profile  and  frequency  bandwidth  obtained 
from  nonlinear,  self-consistent  simulations  are  shown  in  Figs.  8  and  9  respectively.  For  a 
MIG  beam  with  5  %  axial  velocity  spread,  the  FVVJIM  constant-drive  bandwidth  is  6 
with  a  peak  output  power  of  .533  k\V  and  eflicienry  0121.3  %  at  3'l.ri  GHz.  Thes(>  results 
show  that  a  harmonic  gyro-TVVl’  can  deliver  output  power  at  least  an  order  of  magnitude 
higher  than  that  possible  in  state-of-the-art  fundamental  gyrn-TVV'rK. 
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Figure  7:  Oscillation  threshold  lengths  for  T TE[\^  and  TE\\^  gyro-BWO  inodes  as  a 
function  of  beam  current,  is  the  waveguide  radius. 


Table  I:  Second  Harmonic  Gyro-TWT  Design 


Beam  Voltage 

100  kV 

Beam  Current 

25  A 

a  =  ux/t^ll 

1.0 

Magnetic  Field 

6.4  kG 

Center  Frequency 

35  GHz 

Mode 

TE2^ 

Cyclotron  Harmonic 

2nd 

Circuit  Radius 

0.44  cm 

rdr^ 

0.4 

BIB, 

0.98 

Aj-),/r;i| 

5% 

Number  of  Stages 

3 

Sever  Length 

3  cm 

(2) 

Figure  8:  Predicted  spatial  power  profile  of  35  GHz  wave  in  proposed  three-stage  TE^^ 
gyro-TWT  (Table  1)  for  input  power  of  0.5  W  (dashed  curve).  1.9  \V  (solid  curve),  and  4 
W  (dotted  curve). 
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Frequency  ( GHz ) 


Figure  9;  Constant-drive  bandwidth  of  proposed  three-stage  T'F’^i  g.vro-TWT  (Table  I)  for 
axial  velocity  spreads  of  O'X  (dashed  curve),  3%(dotled  curve)  and  bVi.  (solid  curve)  with 
input  power  of  0.21  W,  0.55  VV  and  1.9  \V,  respectively. 

As  seen  from  Fig.  7.  the  start  oscillation  length  for  the  TE\y  gyro-BWO  is  approxi¬ 
mately  three  times  longer  than  for  the  T gyro-BWO.  Therefore,  if  the  TFsi  and  TE\\ 
modes  can  be  eliminated,  a  single-stage  device  can  yield  high  gain  and  yet  remain  stable. 

Such  a  single-stage  proof-of-principle  device  is  being  built  at  L’CL.A  on  the  .4KO  program 
with  the  design  parameters  listed  in  Table  11.  By  slicing  the  cylindrical  waveguide  a.xially 
with  two  cuts  separated  in  azimuth  by  90",  both  the  TE\\  and  7'Fai  modes  will  be  sup¬ 
pressed.  A  similar  procedure  was  successfully  employed  in  a  previous  experiment  [10].  The 
performance  of  the  gyro-TWT  given  in  Table  11  was  also  numerically  simulated.  Compared 
with  the  previous  multi-stage  gyro-TWT  design,  the  single-stage  device  provides  compara¬ 
ble  performance  and  is  less  sensitive  to  the  beam’s  axial  velocity  spread.  For  the  case  of 
8%  axial  velocity  spread,  a  peak  outi)ut  power  of  400  kW.  peak  efficiency  of  20‘7(  and  gain 
of  33  dB  are  predicted  at  15.9  Cllz  and  the  FWHM  constant-drive  bandwidth  is  0  ‘a.  The 
bandwidth  of  the  gyro-TWT  is  shown  in  f  ig.  10. 

Table  II  :  Proof-of-Principle  Experiment  Design 


Beam  Voltage 

100  kV 

Beam  Current 

20  A 

Q  =  I’l  /  t’li 

1.0 

Magnetic  Field 

2.9  kC 

Center  Frequency 

16.1  GHz 

Mode 

TE,: 

Cyclotron  Harmonic 

■2nd 

Circuil  Hadiii.-. 

0.95  cm  i 

r.  /'■„ 

0.4  I 

B/B, 

0.98 

Ar|,/r|, 

8  % 

Circuit  Length 

65  cm 

for  axicil  velocity  spreads  of  0%  (dashed  curve),  3%(dotted  curve)  and  5/c  (solid  curve)  with 
input  power  of  0.21  VV,  0.55  VV  and  1.9  W.  respectively. 


Although  the  interaction  length  has  been  reduced  to  80%  of  the  T gyro-B VVO  thresh¬ 
old  length,  the  gain  of  the  single-stage  device  is  still  appreciable  for  such  a  high  power 
amplifier. 

The  100  kV,  20  A  beam  with  a  =  1  will  be  emitted,  as  shown  in  Fig.  11.  from  the 
curved  cathode  surface  of  the  single  anode  MIG  described  in  Table  Ill.  It  is  predicted  by 
the  EGU.N  code  to  yield  Ar../r-  =  .8%  and  is  being  built  by  .NTHl'.  'Fhe  MIG  beam  will  be 
compressed  adiabatically  into  the  interaction  region  where  a  magnetic  held  with  ±  0.25% 
uniformity  over  1  m  is  provided  by  a  4  kG,  liquid  cooled,  aluminum  foil-wound  solenoid. 
Using  HP’s  high  frequency  structure  simulator,  a  two-mode  axial  phase  velocity  coupler  has 
been  modeled  and  simulated.  The  rf  input  coupler  shown  in  Fig.  12  will  couple  a  rectangular 
waveguide’s  7'E\o  mode  to  a  circular  guide’s  TE2\  mode  with  0  dB  coupling  over  a  20  % 
bandwidth  [11],  The  predicted  bandwidth  of  the  input  coupler  is  shown  in  Fig.  13.  The 
high  power  TE2\  output  wave  from  the  amplifier  will  be  converted  into  tlu'  more  useful 
T E\\  mode  in  the  12  cm  length,  two-period  beat  wave  converter  [11.  12.  13]  shown  m  Fig. 
14  with  a  predicted  0.5  dB  bandwidth  of  10  %.  and  a  conversion  efficiency  of  99.5  %  at  10.1 
GHz.  The  TE21/TE11  mode  converter  design  resulted  from  our  collaboration  with  Dr.  ,1. 
Pretterebner  of  .Stuttgart  University.  The  predicted  performance  of  the  converter  is  shown 
in  Fig.  15. 

KZ)  cruss 


Figure  11:  Trajectories  of  electrons  emitted  from  the  100  k\'.  20  A  single  anode  MKi.  I'he 
electrode  and  magnetic  field  profile  are  also  shown. 
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Table  III  ;  Single  Anode  MIG  Design 


Beam  Voltage 

too  kV 

Beam  Current 

20  A 

Q  =  rj./u|| 

1.0 

Cathode  Loading 

9.7  .A  /  cm’ 

Cathode  Radius 

0.57  cm 

Cathode  Angle 

52“ 

Radius  of  Cathode  Curvature 

4.00  cm 

Emitting  Strip  Length 

0.58  cm 

Cathode  to  Anode  Distance 

4.1  cm 

Magnetic  Compression  Ratio 

3.8 

Ar.'||/e|| 

8 

D.  High  Power  Harmonic  CARMs 


Summary 


The  interaction  strength  of  an  nth  harmonic  TE„i  CARM  amplifier  with  an  axis- 
encircling  electron  beam  can  be  significantly  enhanced  by  reducing  the  magnetic  field 
to  achieve  a  grazing  intersection  of  the  cyclotron -resonance  line  with  the  waveguide 
mode  and  thereby  increasing  the  electron  Larmor  radius  towards  the  wall  radius  where 
the  fields  are  stronger.  A  grazing  intersection  also  considerably  relaxes  the  require¬ 
ments  on  beam  quality.  We  have  performed  preliminary  investigations  of  a  third- 
harmonic  TE31  CARM-TWT  amplifier  whicli  could  be  driven  by  our  3.30  k\  electron 
gun  and  modulator. 

CARMs  [14,  15,  16,  17,  18]  can  generate  very  short  wavelength  emission  much  more 
readily  than  gyrotrons.  Due  to  the  large  axial  velocity  of  the  electrons,  the  output  wave  is 
Doppler-upshifted  by  the  factor  ‘27||  as  in  a  EEL.  For  example,  a  500  k\’,  10  A  C.ARM  in  a 
60  kG  magnetic  field  can  be  expected  to  produce  1  MW  at  500  GHz.  Harmonic  interaction 
results  in  a  further  reduction  of  the  magnetic  field  by  a  factor  of  the  harmonic  number. 
CARMs  can  handle  extremely  high  power  because  the  attenuation  is  weak  since  they  are 
operated  far  above  the  cutoff.  A  similar  voltage  gyrotron  operating  at  500  GHz  would 
require  a  magnetic  field  of  ~  350  kG  and  would  be  limited  to  much  lower  power  levels.  .An¬ 
other  feature  of  C.A.RMs  is  that  they  can  be  extremely  efficient  due  to  the  autoresonance, 
where  the  electrons  stay  in  resoitance  with  the  wave  even  as  they  lose  energy  to  it  due  to 
the  compensation  between  the  two  terms  in  the  resonance  condition.  The  change  is  the 
relativistic  cyclotron  frequency  is  compensated  by  the  change  in  the  Dopplcr  shift .  The 
proof-of-principle  experiment  [19]  that  we  performed  in  collaboration  with  LLNL  produced 
megawatt  power  levels  up  to  270  GHz.  We  are  proposing  in  our  .ARC)  basic  studies  to  care 
fully  provide  the  detailed  physics  basis  for  the  next  generation  of  devices  since  experiments 
to  date  have  not  yet  satisfied  this  need. 

We  have  investigated  the  advantage  of  a  third-harmonic,  multi-section  C'.ARM- I  W  1 
amplifier  driven  by  a  high  voltage,  low  current  axis-encircling  electron  beam  (ilfiO  k\',  2  .-\  1. 
The  axis-encircling  beam  will  be  [)roduccd  by  passing  a  pencil  beam  through  a  gyroresonani 
wiggler.  The  beam  will  emerge  from  the  wiggler  with  a  modest  velocity  ratio  (c^/c.  =  o  1 
of  0.5  and  a  considerable  S;  =:  (1  -  C;/c)‘)”'''‘  of  1.44.  The  resulting  large  value  of  7- 
will  be  used  to  give  a  substantial  enhancentent  of  the  coupling  strength  for  the  harmonic 
interaction.  For  an  n-th  harmonic  interaction,  the  coupling  constant  is  proportional  to  [20] 

Iln  =  Un(9nl^c/^u)]^ 

where  r/,  and  r^.  are  the  Larmor  and  cavity  radii,  respectively,  and  (/„}  is  the  fir,st  zero  of 
J^(x).  Our  ARO-funded  high-harmonic  gyro-devices  with  axis-encircling  beams  that  were 
developed  over  the  last  ten  years  have  been  based  on  this  coupling  constant.  By  using  the 
equations  of  motion  the  equation  for  the  cutoff  frequency  =  </r,ic/r„  ). 

and  the  condition  for  grazing  intersection  of  the  n-th  harmonic  cyclotron  resonance'  line' 
with  the  waveguide  mode  (u-v  =  y^nQ.)  the  above  coupling  constant  can  he  written  as 

Hn  =  I'/nl-Dolx)]*’ 
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where  /Jx  is  tlie  electrons"  transverse  velocity  normalizeti  to  the  speed  of  light.  The  depen¬ 
dence  of  Hn  on  i3_l  for  a  third  harmonic  C'.ARM  TW'l'  is  shown  in  Fig.  Hi(a)  for  several 
values  of  7^.  For  our  proposed  experiment  with  7,  =  1.-4-)  and  Jj_  =  0.3U.  it  can  be  seen 
that  the  interaction  strength  increases  by  a  factor  of  four  compared  to  the  unenlianced  case 
of  7;  =  1.  For  these  parameters,  tlie  interactior.  constant  is  sufficiently  strong  to  yield  a 
high  gciin,  efficient  amplifier. 


Figure  16:  Dependence  of  (a)  tlurd-harmonic  interaction  strength  H3  and  (b)  sixth- 
harmonic  interaction  strength  //e  on  normalized  transverse  velocity  =  v±/c  for  four 
values  of  axial  Lorentz  factor  7j|  =  (1  -  (7|l  =  1.  1.2,  1.44  and  2.0). 

The  marginal  stability  design  procedure,  which  was  recently  developed  and  published 
by  our  group  in  collaboration  with  Drs.  .^.T.  Lin  and  K.R.  Chu  for  creating  stable,  high 
performance  gyro-TVVTs  [4]  can  be  used  to  design  the  C.ARM-TVVT  circuit.  In  this  design 
procedure,  the  amplifier  is  comprised  of  several  interaction  sections  separated  by  attenuating 
severs  and  each  with  a  length  shorter  than  the  start-oscillation  length  for  the  most  prevalent 
oscillation  mode,  which  is  usually  a  gyro-BWO  interaction.  The  circuit  will  be  designed  for 
a  350  kV,  2  A  electron  beam  with  an  o  of  0.5  by  finding  the  start-oscillation  length  for  all 
the  possible  modes  and  then  setting  the  length  for  the  interaction  sections  to  be  ~  109(  less 
than  this.  The  desired  gain  for  the  amplifier  determines  the  number  of  interaction  sections 
needed.  The  cylindrical  waveguide  circuit  will  be  sliced  axially  with  three  cuts  separated  by 
120“^  in  azimuth  to  suppress  all  modes  except  tho.se  with  a  threefold  azimuthal  symmetry. 
Figure  16(b)  shows  the  coupling  constant  for  the  strongest  surviving  ntode.  the  TEei  mode. 
Since  the  value  of  He  is  smaller  than  for  our  design  mode  by  a  factor  of  ten.  we  can  expect 
significant  gain  for  the  relatively  long  interaction  length  available. 

A  CARM’s  desirable  features  of  high  efficiency,  high  power,  and  large  output  frequency 
to  rest-mass  gyrofrequency  ratio  can  be  attained  only  if  the  electron  beam's  axial  velocity- 
spread  can  be  kept  low  ( <  2%),  which  is  difficult  to  achieve  for  a  high  current  beam 
with  the  requisite  transverse  velocity.  VVe  have  investigated  gyroresonant  wigglers  [21]  for 
the  production  of  low  velocity-spread  helical  beams.  Using  a  PIC  simulation  code,  we  found 
that  a  near-optimized  gyroresonant  wiggler  designed  for  a  final  o  of  O.G  resulted  in  an  axial 
velocity  ratio  of  6.5%  for  a  150  .\  beam  as  shown  in  Fig.  17  and  4%  for  20  A.  Emitlatice 
growth  is  due  to  the  space-charge  of  the  non-axisymmetric  helical-shaped  beam.  These 
numbers  were  found  using  a  self-consistent  beam  simulation  code  which  in^'ludes  space- 


Figure  17;  Spatial  growili  of  axial  velocity-spread  due  to  gyroresonant  wiggler  at  z  =  0  cm 
to  z  =  40  cm  (250  kV;  150A;  beam  radius  =  1  cm;  wall  radius  =  2  cm:  wiggler  period  =  8 
cm:  relative  wiggler  amplitude.  B^.fBo  =  0.02). 


charge  and  which  wa.s  provided  to  us  by  our  collaborator  Dr.  M.  Caplan  of  LL.N'L.  This 
code  will  be  used  to  design  the  wiggler  for  the  proposed  experiment.  We  expect  to  be  able 
to  achieve  an  axial  velocity  spread  better  than  2%  for  our  2  A  beam.  In  addition,  since  the 
initial  experiment  will  be  performed  with  a  grazing  intersection,  the  beam  velocity-spread 
requirements  are  considerably  relaxed. 

The  parameters  of  a  our  proof-of-principle  experiment  are  given  in  Table  IV.  It  will 
utilize  equipment  available  in  our  laboratory.  Our  high  voltage  SLAC  5045  Klystron  gun 
driven  by  our  SLAC-style  400/800  k\’  modulator  will  create  the  initial  pencil  beam.  The 
experiment  w'ill  be  performed  in  our  3  kG.  2  m  long.  20  cm  bore  copper  solenoid.  The 
vacuum  system  is  currently  operating  in  the  10“**  Torr  range. 

Table  IV.  Parameters  of  33  GHz  third-harmonic  sliced  TEaj  CARM-TWT 


Voltage 

350  kV 

Current 

2  A 

Q  =  rj./n|| 

0.5 

Power 

140  kW 

Efficiency 

2091 

Magnetic  Field 

3.3  kG 

Cluiding  Center  Radius,  7^, 

0 

Wall  Radius, 

0.88  cm 

Larmor  Radius,  r/. 

0.32  cm 

Wiggler  Period 

4.5  cm 

Initial  Beam  Radius 

O.C  mm 
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E.  High  Performance  Gyro-Klystron 


Summary 


The  objective  of  this  work  is  to  investigate  the  basic  physics  issues  associated  with  the 
design  and  operation  of  a  high  performance  millimeter- wave  gyroklystron  amplifier. 
Our  approach  is  to  begin  by  improving  the  performance  of  V’arian’s  pioneering  28 
GHz  gyroklystron  by  employing  our  physics  design.  A  state-of-the-art  self-consistent 
nonlinear  simulation  code  to  design  the  amplifier  was  provided  to  us  by  our  colleague 
Dr.  K.R.  Chu.  Using  the  simulation  code,  a  2.50  kW,  33  GHz  gyro-klystron  with  an 
efficiency  of  39%  and  saturated  gain  of  52  dB  has  been  designed  for  an  80  kV,  8  A 
electron  beam  with  o  =  1.5  and  an  tixial  velocity  spread  of  7%.  A  superconducting 
solenoid  appropriate  for  the  experiment  was  donated  by  Hughes  EDD  and  several  28 
GHz  gyrotron  MIG  electron  guns  were  donated  by  Varian  Associates. 

Our  initial  work  on  the  ARO  program  has  been  on  the  development  of  high  perfor¬ 
mance,  high  power  millimeter- wave  fundamental-frequency  gyro-klystrons.  However,  it  is 
envisioned  that  our  research  will  soon  lead  to  the  development  of  high-harmonic  gyro- 
klystrons  using  high  power  MIG  electron  beams.  Whereas  our  previous  fifth-harmonic  [22] 
gyro-klystron  used  a  more  moderate  power,  gyroresonantly  rf-accelerated  axis-encircling 
electron  beam  and  generated  0.5  kW.  our  future  harmonic  gyro-klystrons  will  be  at  the  100 
kW  output  power  level.  Harmonic  operation  allows  the  magnetic  field  to  be  reduced  by  the 
harmonic  number.  A  third-harmonic,  100  GHz  gyro-klystron  would  require  a  magnetic  field 
of  only  14  kG.  Because  a  gyro-klystron’s  cavities  are  fairly  short  in  order  to  avoid  oscilla¬ 
tion,  operation  of  a  harmonic  gyro-klystron  is  much  more  straightforward  than  operation 
of  a  harmonic  gyrotron.  We  therefore  are  confident  that  200  GHz  gyroklystrons  could  be 
readily  developed  and  have  this  as  our  long  term  aim. 

A  self-consistent  computer  [23]  code  written  by  Professor  K.R.  Chu  and  used  both  for 
the  University  of  Maryland’s  high  power  gyroklystron  and  our  large-orbit  gyro-klystron 
was  employed  to  optimize  parameters  for  our  first  gyro-klystron,  which  is  intended  to  be  a 
significant  extension  of  Varian  .Associates’  pioneering  28  GHz  gyroklystron  [24],  The  code 
allows  for  beam  loading  of  the  cavities,  beam  velocity  spread  and  radial  guiding  center 
spread.  Additionally,  the  user  may  choose  the  analysis  for  each  cavity  to  be  either  linear 
or  nonlinear.  The  code  was  used  successfully  in  predicting  the  results  of  the  University 
of  Maryland’s  gyroklystron  and  our  ARO-funded  fifth  harmonic,  large-orbit  gyroklystron. 
The  code  was  also  tested  against  Varian’s  earlier  results  with  good  agreement.  Varian 
designed  a  two-cavity  28  GHz  gyroklystron  for  ECH  applications  with  a  circular  TEqi  input 
cavity,  a  circular  TE02  output  cavity  and  an  a  =  rx/uj|  =  2.  Severe  oscillation  problems 
were  encountered,  but  after  some  redesign,  an  output  power  of  50  kW  was  obtained.  The 
efficiency  remained  under  10%..  Figure  18  shows  the  nonlinear  gain  predicted  by  the  above 
mentioned  nonlinear  code  for  the  original  Varian  experimental  device  for  varying  the  output 
cavity  Q,  which  allows  us  to  determine  the  peak  efficiency  by  searching  for  the  electric  field 
amplitude  which  gives  optimum  energy  extraction  Figure  19  shows  the  dependence  of 
nonlinear  gain  on  the  beam’s  axial  velocity  spread.  It  is  evident  that  the  efficiency  for  this 
design  is  highly  dependent  on  the  velocity  spread.  This  explains  why  tlie  Varian  device 
displayed  a  rather  poor  efficiency. 
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Figure  18:  Gain  characteristics  of  Varian  two-cavity  gyro-klystron  for  several  values  of 
output  cavity  Q. 
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Figure  19:  Gain  characteristics  of  V’arian  two-cavity  gyro-klystron  for  several  values  of  axial 
velocity  spread. 


The  purpose  of  our  new  experiiiieiu  is  to  improve  V'arian's  earlier  gyroklyslron  and 
produce  a  high  power,  high  performance  amplifier.  An  additional  buncher  cavity  has  been 
added  to  increase  gain  and  the  design  was  completed  with  a  lower  a  for  better  stability. 
The  drift  tubes  will  be  loaded  with  a  thin  ceramic  tube  with  a  lossy  inner  coating  cis  used 
for  the  absorbing  sever  in  a  recent  successful  gyro-TVVT  experiment  [25].  The  frequency 
has  been  increased  to  33.2  GHz  and  the  magnetic  field  was  optimized  for  the  new  frequency. 

The  operating  parameters  of  the  new  design  along  with  simulation  results  can  be  seen  in 
Table  V.  For  an  ideal  beam,  the  code  predicts  an  optimum  output  power  of  300  kW,  51  dB 
saturated  gain  and  46%  efficiency.  Using  a  more  realistic  axial  velocity  spread  of  7%,  the 
predicted  results  are  an  output  power  of  250  kVV  with  an  efficiency  of  39%  and  a  saturated 
gain  of  49  dB.  Figure  20  shows  the  predicted  linear  and  nonlinear  gain  characteristics  for 
several  values  of  velocity  spread.  It  can  be  seen  that  the  device  will  saturate  for  an  input 
power  of  approximately  6  W. 


Table  V.  Parameters  of  High  Performance,  Fundamental  Three-Cavity 

Gyro-Klystron  Amplifier 


Voltage 
Current 
Frequency 
Magnetic  Field 
a  =  rj./u|| 

Mode'' (#1,  #2.  #3) 
Q  (#1,#2,  #3) 
f^cavity  (#1.  #2,  #3) 

Tdrift-tub«-  (#1,  #2) 

Circuit  Length 


80  kV 
8  A 

33.2  GHz 
13  kG 
1.5 
7% 

0.48 

TEoi,  TEoi,  TEo2 
300, 300, 550 
2A,  2A,  2A 
3A,  3A 
10.8  cm 


Simulation  Results: 

Power  250  kW 

Efficiency  39% 

Gain  (Saturated)  52  dB 


The  electron  beam  will  be  produced  by  a  28  GHz  gyrotron  MIG  donated  by  Varian 
Associates.  A  32.9-33.5  GHz  tunable.  50  kW  magnetron  will  initially  provide  the  input 
signal.  Eventually,  this  will  be  replaced  by  a  10  W  helix  TWT  Ka-band  amplifier.  A  seven- 
coil  superconducting  solenoid  with  a  12  cm  uniform  region  donated  by  Hughes  Aircraft  Co. 
will  be  used  to  generate  the  12.4  kG  field  necessary  in  the  interaction  region.  All  three 
cavities  are  2A  long  and  are  connected  by  3A  long  drift  tubes.  The  driver  and  buncher 
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Figure  21:  CAD  layout  of  gyroklystron  circuit  within  its  vacuum  chamber  and  supercon¬ 
ducting  solenoid  dewar. 
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cavity  will  operate  in  the  TEqi  mode  and  the  output  cavity  operates  in  the  TE02  mode  so 
as  to  allow  cw  operation.  The  circui*  layout  has  been  completed  using  a  CAD  program. 
The  full  drawing  is  shown  in  Fig.  21. 

The  diffraction  coupled  output  cavity  was  designed  using  a  computer  code  written  by 
Professor  K.R.  Chu  for  the  UCLA  ATRl  program  [26].  Given  the  physical  dimensions  of 
the  cavity  and  the  TE^n/  mode  in  which  the  cavity  is  to  operate  in,  the  code  predicts  the 
field  distribution  in  the  cavity,  the  cavity’s  Q,  and  the  resonant  frequency.  The  code  works 
for  any  arbitrarily  shaped  cylindrical  cavity  of  slowly  varying  cross-section.  The  cavity 
must  be  closed  at  one  end  (the  cross-section  must  be  smaller  than  the  cut-off  cross-section), 
and  can  be  either  open  or  closed  at  the  other  end.  The  code  can  also  account  for  wall 
loss  in  the  cavity.  As  in  the  Varian  Associates  experiment,  it  was  decided  to  design  the 
cavity  for  TE021  operation  with  a  length  of  2A.  The  cavity  was  designed  to  have  a  Q  of  550 
and  an  operating  frequency  of  33.2  GHz  as  specified  by  the  gyroklystron  simulation  code. 
The  profile  of  the  output  cavity  is  shown  in  Fig.  22(a)  and  Fig.  22(b)  shows  the  resultant 
electromagnetic  profile  of  the  cavity’s  TE021  mode. 
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Figure  22:  (a)  Dependence  of  inner  wall  radius  on  axial  distance  for  gyroklystron  diffraction- 
coupled  output  cavity  and  (b)  axial  profile  of  gyroklystron  output  cavity's  TE021  mode 
including  emitted  wave. 

To  prevent  oscillation  from  occurring  in  the  beam  compression  region,  lossy  MgO-SiC 
rings  have  been  ordered  and  received.  To  keep  the  drift-tubes  between  the  gyroklystron 
cavities  from  oscillating,  we  will  use  a  technique  successfully  demonstrated  in  a  recent  gyro- 
TWT  experiment  [25].  An  excellent  attenuator  can  be  achieved  with  a  thin  dielectric  tube 
coated  with  graphite  on  the  inside  wall.  In  addition,  the  input  and  gain  cavities  must  also 
be  loaded  with  lossy  rings  in  order  to  achieve  the  correct  value  of  Q  in  the  normally  high-Q 
TEoi  cavities.  Two  options  are  available  for  this.  Because  carbon  impregnated  aluminum 
silicate  rings  [27]  have  the  advantage  that  the  turn-around  time  can  be  relatively  quick,  we 
have  devoted  considerable  effort  to  develop  this  technique.  The  other  option  is  to  again  use 
lossy  carbon  coated  dielectric  tubes  for  this  application  too. 

The  magnetic  field  required  for  the  MIG  electron  gun  will  be  produced  by  copper 
solenoids  external  to  the  superconductor’s  dewar.  The  profile  has  been  designed  to  con¬ 
form  to  the  magnetic  field  profile  of  the  Varian  gyrolron.  Our  magnetic  field  profile  is 


Figure  23:  Gyro-klyslron  s  magnetic  field  profile  on  axis  with  and  without  external  gun 
coils  activated.  The  position.^  of  the  superconducting  and  external  coils  are  shown  as  well 
as  the  cathode  location. 

shown  in  Fig.  23.  The  solenoids  are  currently  being  mounted  to  the  superconducting  de- 
war.  The  vacuum  chamber  and  pump.s  have  been  assembled  and  produce  a  vacuum  in  the 
10"  range.  Also,  we  have  constructed  a  modulator  to  drive  the  50  kW,  33.2  GHz  tun¬ 
able  magnetron  that  will  serve  as  the  initial  input  source  for  the  gyroklystron  amplifier. 
A  high  voltage  power  supply,  thyratron  and  pulse  transformer  were  assembled  and  a  40  f) 
distributed  transmission  line  was  built  and  tested.  The  magnetron  was  driven  to  its  full  50 
k\\  output  power.  Directional  couplers  and  attenuators  are  available  to  reduce  the  power 
to  the  required  drive  level. 

Penultimate  cavity  detuning  [2S.  29]  has  also  been  investigated.  Similar  to  the  linear 
klvstron,  a  gyroklystron  s  efficiency  may  be  enhanced  by  properly  misluning  the  bunching 
cavities  so  as  to  increase  the  electron  bundling.  In  the  gyroklystron,  it  is  desirable  to  have 
the  leading  phase  electrons  accelerated  and  lagging  phase  electrons  decelerated.  This  can 
be  accomplished  by  detuning  the  cavity  frequency  downward.  It  can  be  seen  in  Fig.  24  that 
for  zero  velocity  spread,  significant  enhancement  of  the  efficiency  is  predicted,  but  at  the 
expense  of  gain.  This  loss  in  gain  could  be  compensated  by  the  addition  of  more  bunching 
cavities.  However,  when  a  more  realistic  velocity  spread  of  was  included,  simulation 
showed  that  the  increase  in  efficiency  was  negligible  and  the  gain  was  still  greatly  reduced, 
for  this  reason,  detuning  of  the  penultimate  cavity  in  our  three-caviiy  gyro-klystron  is  not 
planned. 
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Figure  24:  Efficiency  and  saturated  gain  of  UCLA  gyro-klystron  as  a  function  of  penultimate 
cavity  detuning  for  two  value  of  axial  velocity  spread  =  0%  and  7%). 
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Marginal  stability  design  criterion  for  gyro-TVVTs  and  comparison  of 
fundamental  with  second  harmonic  operation 

A  T.  LINt.  K.  R.  CHUt:.  C.  C  LlNt.  C.  S.  KOU??.,. 

D  B.  McDERMOTT§  and  N.  C  LUHMANN  Jr§ 

Stability  properties  of  both  the  fundamental  and  second  harmonic  gyrotron 
travelling  wave  amplifier  (gyro-TWT)  arc  examined  with  multi-mode  particle 
simulations.  The  second  harmonic  cyclotron  interaction  with  an  axis-encircling 
electron  beam  is  found  to  be  more  stable  to  oscillations  and  can  yield  significantly 
greater  power  than  the  fundamental  harmonic  gyro-TWT  A  multiple  stage 
interaction  structure  based  on  a  marginal  stability  criterion  is  proposed  and 
illustrated  with  examples  of  a  !28kW  fundamental  gyro-TWT  and  a  532  kW 
second  harmonic  gyro-TW'T  Stable  amplification  at  much  higher  power  levels  is  in 
principle  possible. 


1.  Introduction 

The  gyrotron  travelling  wave  amplifier  has  superior  power  handling  capability 
compared  to  the  conventional  TWT  because  the  gyro-TWT  employs  a  fast  wave 
interaction  and  a  smooth  wall  geometry.  In  the  millimetre-wave  band,  the  demon¬ 
strated  peak  power  (~25kW)  (Barnett  et  al.  1979,  1989.  Granatstem  ei  at.  1983. 
1984.  Chu  el  al.  1990  a)  and  bandwidth  (Chu  ei  al.  1990  a)  {~  10%)  of  gyro-TWTs 
are  already  comparable  to  those  of  state-of-the-art  TWTs  (see,  for  example, 
Butterworth  and  Wallace  1987).  However,  in  these  demonstrations  of  the  fundamen¬ 
tal  harmonic  gyro-TWT,  a  relatively  large  magnetic  field  (~  I2kG  for  the  Ka-band) 
was  required.  Also,  the  power  levels  achieved  are  still  too  low  fully  to  realize  the 
power  handling  capability  of  the  gyro-TWT. 

The  strong  beam-wave  interaction  at  the  fundamental  cyclotron  harmonic  has 
been  the  cause  of  unwanted  oscillations  which  limit  the  output  power.  There  arc 
various  ways  to  increase  the  oscillation  threshold,  such  as  (i)  lowering  the  beam 
a  {  =  vjv^)  to  reduce  the  interaction  strength,  (ii)  increasing  the  beam  voltage  to 
move  away  from  the  cut-ofT  frequency  of  the  operating  mode,  and  (iii)  employing  a 
severed  interaction  waveguide.  While  all  of  these  methods  have  proved  cfTective  in 
improving  the  stability  and  power,  additional  measures  are  required  if  still  higher 
power  is  desired. 

In  this  study,  wc  investigate  methods  of  stably  operating  the  fundamental  and 
harmonic  gyro-TWT  as  a  high  power  amplifier.  Harmonic  interaction  is  much 
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weaker,  but  it  ofTers  the  important  advantage  of  a  reduced  magnetic  field  require¬ 
ment.  To  compensate  for  the  weaker  interaction,  one  can  introduce  an  azimuthally 
corrugated  interaction  structure  (Lau  and  Barnett  1982,  Chu  and  Dialetis  1985, 
Namkung  1984.  Destlcr  1983)  which  ripples  the  RF  field  lines  in  a  way  to  strengthen 
a  specific  harmonic  field  component.  Alternatively,  one  can  employ  a  large  orbit 
concentric  electron  beam  (hundreds  of  kilovolts)  (Destler  1983,  McDermott  et  al. 
1983.  Kou  ei  al.  1990)  with  the  electrons  located  in  the  vicinity  of  the  peak  field  of  an 
azimuthal  mode.  However,  the  brute-force  remedy  of  increasing  the  electron  beam 
current  (as  considered  in  this  study)  may  be  the  most  suitable  for  high  power 
operation.  It  ofTers  the  dual  merits  of  moderate  voltage  and  smooth  waveguide 
structure. 

Our  study  is  concerned  with  power  levels  well  beyond  the  state  of  the  art.  Self- 
oscillations  are  inevitability  present  in  such  a  regime.  It  is  thus  essential  that  a  multi¬ 
mode  analysis  be  employed  to  examine  the  physical  behaviour  of  these  modes.  More 
importantly,  means  of  stabilizing  the  unwanted  oscillations  must  also  be  investigated 
so  as  to  push  upward  the  stability  boundary  and  hence  the  power  output. 

The  basic  tools  of  the  subsequent  studies  consist  of  an  analytical  linear  theory 
(Chu  and  Lin  1988),  a  single-mode  particle  tracing  code  (Kou  et  al.  1991),  and  a 
multi-mode  particle  simulation  code  (Lin  and  Lin  1989),  which  complement  and 
validate  one  another.  The  first  two  tools  allow  us  to  determine  the  parameters  for 
optimum  single-mode  operation.  The  resultant  stability  properties  are  then  exam¬ 
ined  by  the  multimode  panicle  simulation  code. 

In  §2,  the  simultaneous  evolution  of  various  oscillating  modes  are  analysed  in 
multi-mode  simulation,  which  leads  to  the  conclusion  that  the  second  harmonic 
gyro-TWT  is  generally  more  stable  than  the  fundamental  harmonic  gyro-TWT. 
Based  on  these  stability  considerations,  a  procedure  is  developed  in  §3  for  the  design 
of  a  marginally  stable  gyro-TWT.  It  is  found  that,  under  the  marginally  stable 
operating  condition,  the  second  harmonic  gyro-TWT  produces  significantly  greater 
power  than  the  fundamental  harmonic  gyro-TWT  with  comparable  gain,  band¬ 
width.  and  efficiency. 


2.  Self-oscillation  in  fundamental  and  second  harmonic  gyro-TWTs 

Start-oscillation  conditions  for  the  forward  (Chu  and  Lin  1988,  Lau  et  al.  1981, 
Davies  1989)  and  backward  (Park  et  al.  1984)  waves  in  the  gyro-TWT  and  cyclotron 
autoresonance  maser  amplifiers  have  been  extensively  studied  in  the  literature. 
Theoretical  calculations  and  experimental  evidence  all  indicate  that  oscillation  of 
one  kind  or  another  is  the  norm  rather  than  exception  in  a  high  power  gyro-TWT.  In 
this  section,  we  investigate  the  spontaneous  growth  of  these  oscillations  and  their 
saturated  behaviour  through  particle  simulations.  To  determine  the  appropriate 
parameters  to  use  for  the  simulation,  we  have  employed  the  linear  theory  (Chu  and 
Lin  1988)  and  the  single-mode  particle  tracing  code  (Kou  et  al.  1991)  to  arrive  at  the 
two  cases  shown  in  the  Table.  The  electron  beam  in  Case  a  is  consistent  with  a 
magnetron  injection  gun  (Baird  and  Lawson  1986),  while  Case  b  corresponds  to  a 
cusp  electron  gun  (Scheitrum  et  al.  1989).  Each  case  represents  a  compromise 
between  stability,  gain,  power,  bandwidth,  and  efficiency.  Here,  is  the  magnetic 
field  which  results  in  a  grazing  intersection  between  the  operating  mode  and  the 
operating  cyclotron  harmonic  resonance  line.  In  Fig.  I,  the  operating  mode  is  shown 
by  the  (near  grazing)  intersection  of  the  two  unbroken  curves.  Intersections  that 
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occur  for  negative  k,  arc  potential  sources  of  oscillations.  For  Case  b  however,  only 
the  5  =  m  modes  arc  unstable  because  the  beam  is  assumed  to  be  axis-encirelmg. 
Furthermore,  all  the  possible  oscillating  modes  have  a  harmonic  number  equal  to  or 
greater  than  the  operating  harmonic  (s  =  2). 

To  emphasize  and  enhance  the  self-oscillations,  we  first  perform  the  simulations 
in  a  relatively  long  (15  cm,  unsevered)  interaction  waveguide  in  the  absence  of  an 
input  signal.  Figure  2  shows  the  temporal  behaviour  of  the  TE^,,  field  amplitudes 
(superscript  indicates  the  harmonic  number  of  the  interaction)  from  noise  to 
saturation.  There  are  three  spontaneously  growing  modes  for  the  case  of  fundamen¬ 
tal  harmonic  interaction  (Case  a)  with  the  operating  mode  being  the  dominant  one 
at  saturation.  In  contrast,  only  one  spontaneously  growing  mode  is  present  for  the 
second  harmonic  interaction  (Case  b)  despite  a  five-fold  increase  in  beam  current.  It 
is  also  significant  that  the  operating  mode  is  stable  with  respect  to  the  absolute 
instability  in  Case  b.  Comparison  between  Case  a  and  Case  b  supports  our  earlier 
assertion  that  the  harmonic  interaction  is  a  viable  scheme  for  high  power  generation, 
This  point  will  be  further  reinforced  by  efficiency  and  band^sidth  considerations  in 
the  following  section. 

To  quantify  the  effects  of  these  oscillations.  Figs  3  and  4  display,  respectively,  the 
electron  momentum  and  energy  distributions  corresponding  to  the  saturated  state  of 
Fig.  2.  The  beam  quality  is  so  degraded  that  no  significant  amplification  can  be 
expected  to  take  place.  As  an  illustration.  Fig.  5(a)  shows  the  temporal  growth  of  a 
TE21  drive  wave  amplified  by  the  spoiled  beam  described  by  Figs  3(/))  and  4(h). 
Figure  5(b)  shows  the  temporal  growth  of  the  same  drive  wave  due  to  a  good  quality 
beam  (see  the  Table,  Case  b)  with  the  oscillating  modes  artificially  suppressed  (i.e. 
single-mode  model).  Comparison  of  Figs  5(o)  and  5(b)  dramatically  demonstrates 
the  deficiency*of  a  single-mode  treatment  of  high  power  gyro-TWTs. 

3.  Marginal  stability  design  procedure  for  high  power  generation 

The  dominance  of  self-oscillations  in  the  cases  examined  in  the  previous  section 
suggests  that  measures  must  be  taken  to  stabilize  these  oscillations  in  order  for  the 
gyro-TWT  to  function  properly  as  an  amplifier.  There  exists  a  tuning  range  for  the 
magnetic  field  of  at  most  a  few  percent  for  normal  gyro-TWT  operation.  While  this 
is  a  sufficient  tuning  range  to  stabilize  the  self-oscillation  of  the  operating  mode,  it 
offers  no  remedy  for  other  oscillating  modes.  Reducing  the  beam  current  can 


Case  a 

Case  b 

Waveguide  radius  (r,) 

0-26  cm 

0  44  cm 

Operating  mode 

TE., 

TE,, 

Operating  cyclotron  harmonic  (s) 

1 

n 

Beam  voltage  (V,,) 

100  kV 

100  kV 

Beam  current  (1^) 

5  A 

25  A  • 

Guiding  centre  position  (r,) 

0  35  r. 

0 

Electron  velocity  ratio  ai(v^otv,„) 

10 

10 

Axial  velocity  spread  (Av,/v,(,) 

5% 

Magnetic  field  (Bq/B,) 

09K 

0-98 

Parameters  of  the  fundamental  harmonic  (Case  a)  and  sexond  harmonic  gyro-TWT  (Case  h) 
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eventually  stabilize  all  oscillations,  but  this  is  normally  associated  with  reduced 
output  power.  The  selective  suppression  of  undesirable  modes  by  special  circuit 
design  can  be  effective  in  some  cases,  but  does  not  have  universal  applicability 
Because  of  the  sensitivity  of  the  oscillation  threshold  to  the  interaction  length 
(Davies  1989),  reducing  the  interaction  length  appears  to  offer  a  practical  and 
effective  means  for  oscillation  stabilization.  The  resultant  loss  in  gain  can  be 
recovered  with  a  multi-stage  interaction  structure  Indeed,  a  recent  severed  gyro- 
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Figure  2.  Field  amplitudes  of  self-oscillating  modes  versus  time  in  an  unsoscred  inicrucimn 
waveguide  of  15cm  in  length  for  («)  Case  a  and  (/>)  Case  h  (supcrscnpi  mdicaics 
cyclotron  harmonic  number). 
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Figure  3.  Distribution  of  transverse  electron  momentum  at  saturation  for  (o)  Case  a  and  (6) 

Case  b. 


TWT  experiment  has  convincingly  demonstrated  (Chu  ei  al.  1990  b)  the  effectiveness 
of  such  a  scheme.  Below  we  explore  a  design  procedure  based  on  this  approach 
which  can  be  systematically  employed  for  high  power  gyro-TWT  designs. 

The  design  procedure  is  based  on  a  marginal  stability  criterion,  namely  the 
conditions  of  each  single  stage  are  pushed  to  the  point  of  marginal  stability  at  the 
desired  beam  power  and  multiple  stages  (separated  by  attenuators)  are  employed  to 
provide  the  desired  gain. 

To  implement  the  marginal  stability  design,  we  first  define  a  critical  waveguide 
length  (Lj)  as  the  length  below  which  all  the  modes  arc  stable  and  above  which  one 
or  more  modes  will  be  spontaneously  oscillating.  This  length  can  readily  be  obtained 
by  iterative  simulation  runs  with  all  the  likely  oscillating  modes  included.  The 
critical  length  thus  obtained  is  observed  in  Fig.  6  to  be  a  sensitive  function  of  the 
beam  current.  At  beam  currents  of  interest  to  high  power  generation,  Z..  becomes 
too  short  to  achieve  significant  gain  in  a  single  section.  Hence,  multiple  stages  are 
required.  The  simulation  results  provide  a  useful  design  guide.  Upon  fixing  the  beam 
current  from  power  considerations,  the  marginal  waveguide  length  can  be  found 
from  Fig,  6,  which  is  then  taken  to  be  the  length  of  each  section  of  the  multi-stage 
gyro-TWT.  Since  each  stage  is  stable  to  oscillations,  the  single-mode  particle  tracing 
code  (Kou  ei  al.  1991)  can  be  employed  for  the  subsequen,  optimization  of 
bandwidth,  gain,  and  efficiency. 

Figures  7  and  8  present  examples  of  multi-stage  gyro-TWT  designs  using  the 
parameters  of  the  Tabic.  For  Case  a,  the  interaction  length  must  be  reduced  to 
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4  5  cm  completely  to  avoid  oscillation  of  the  various  modes.  To  obtain  a  substantial 
gain,  a  three-stage  system  is  required.  Simulation  results  of  the  fundamental  gyro- 
TWT  based  on  the  single  mode  particle  tracing  code  arc  shown  in  l  ig  7  Employing 
a  5  A.  lOOkV  electron  beam  with  an  axial  velocity  spread  of  5"o.  the  fundamental 
harmonic  gyro-TWT  is  predicted  to  achieve  a  maximum  saturated  gain  of  50 dB, 
saturated  output  power  of  128kW  (25  6%  efiiciency).  and  a  .7  2 GHz  bandwidth 
(9-3%)  On  the  other  hand,  an  interaction  length  of  7-95 cm  sufTices  to  exclude  all 
oscillations  for  Case  b.  A  three-stage  system  is  again  used  to  provide  ample  gain. 
Simulation  results  of  the  second  harmonic  gyro-TWT  are  presented  in  Fig.  8  Using 
a  25  A,  lOOkV  electron  beam  with  an  axial  velocity  speed  of  5'’o.  the  second 
harmonic  gyro-TWT  is  predicted  to  achieve  a  maximum  saturated  gam  of  54 dB.  a 
saturated  output  power  of  532kW  (21-3%  efiiciency).  and  a  2  OGHz  bandwidth 
(8-7%). 

4.  Summary  and  discussion 

This  investigation  begins  with  a  multi-mode,  time-dependent  simulation  study  of 
the  stability  properties  of  the  gyro-TWT.  For  high  beam  power,  it  is  shown  that  self- 
oscillations  arising  from  the  interactions  of  various  waveguide  modes  and  dilTerent 
cyclotron  harmonics  can  grow  to  significant  saturation  le\els  such  that  they 
completely  spoil  the  amplification  process.  To  address  the  critical  issue  of  stability, 
we  propose  a  marginal  stability  design  procedure  which  has  been  motivated  by  two 
recent  experimental  developments.  The  demonstration  of  a  severed  gyro-TWT  (Chu 
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Figure  4  Distribution  of  electron  energy  at  saturation  for  [a)  Case  a  and  (/>)  Case  h 
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Figure  5.  Temporal  growth  of  TEji  drive  wave  for  Case  b  due  to  (o)  spoiled  beam  of  Figs 

3(6)  and  4(6)  and  (6)  good  quality  beam. 


et  al.  1990  b)  proves  the  feasibility  of  a  multi-stage  interaction  structure  which  is  an 
essential  element  of  the  design  scheme.  The  demonstration  of  the  suppression  of 
oscillation  by  a  drive  wave  (Barnett  ei  al.  1989,  Chu  et  al.  1990  a,  1991)  justifies  the 
proposed  borderline  design  concept,  since  weak  oscillations,  if  present,  can  be 
cfTeclively  suppressed  by  the  drive.  In  many  applications,  where  occasional  oscil¬ 
lation  is  unacceptable,  an  unconditionally  stable  gyro-TWT  can  be  obtained  by 
choosing  the  section  lengths  to  be  less  than  L,  by  a  safety  margin  of  ~5%. 

The  stability  analysis  and  design  examples  all  indicate  that  a  gyro-TWT 
operating  at  the  second  cyclotron  harmonic  has  significantly  greater  power  capa¬ 
bility  than  the  fundamental  harmonic  version.  As  one  keeps  increasing  the  output 
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power  through  the  marginal  stability  design  procedure,  an  ultimate  limit  will  be 
imposed  by  factors  not  accounted  for  in  the  present  model.  DC  space  charge 
depression  of  the  electron  beam  is  a  well-known  limiting  factor  m  beam  formation. 
AC  space-charge  effects  (Charbit  et  al.  1981,  Liu  and  Yang  1981.  Chen  and  Chu 
1986.  Chu  and  Lyu  1986,  Antonsen  ct  al.  1986.  Latham  1990)  can  modify  the 
operating  mode  (growth  rate,  mode  structure,  etc.)  or  introduce  entirely  new 
unstable  modes.  The  ultimate  capacity  of  the  gyro-TWT  will  not  be  known  until  we 
have  a  belter  understanding  of  these  various  high  current  effects. 

Finally,  several  caveats  on  these  high  power  gyro-TWTs  should  be  discussed.  In 
the  analysis,  the  ends  of  each  section  were  assumed  to  be  perfect  matches,  not  only 
at  the  operating  frequency,  but  also  for  the  various  harmonic  interactions  This  is 
especially  difficult  to  achieve  for  the  input  coupler.  The  gyro-BWO  stability  length 
can  be  expected  to  decrease  for  finite  reflections.  Fortunately,  the  relevant  parameter 
for  oscillation  is  the  product  of  the  two  reflection  coefTicienls  irul  it  is  fairlv 
straightforward  to  build  a  broadband  attenuator,  allhon^h  somewhat 

complicied  by  the  fact  that  the  most  serious  competing  >  \hc  fii.M  atid 
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Figure  6.  Dependence  of  start  oscillation  length  on  beam  current  for  («)  Case  a  and  (6i 

Case  b. 
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Frequency  (GHz) 

Figure  7.  Performance  of  Case  a  using  a  three-stage  design:  (o)  Spatial  power  profile  of 
34  6  GHz  wave  driven  into  saturation  and  (6)  saturated  output  power  and  efficiency 
versus  frequency. 

second  harmonic  amplifiers  threatens  to  oscillate  near  its  cut-olT.  A  further  concern 
is  the  high  power  capability  of  the  attenuating  sever.  The  attenuator  used  in  the 
recent  severed  gyro-TWT  experiment  (Chu  et  al.  1990  b)  was  a  tapered  ceramic  tube 
with  a  lossy  inner  coating.  An  improved  version  may  have  to  be  designed  before 
these  devices  are  operated  with  high  average  power. 
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Figure  8.  Performance  of  Case  b  using  a  three-stage  design:  (o)  Spatial  power  profile  of 
34  GHz  wave  driven  info  saturation  and  (b)  saturated  output  power  and  efficiency 
versus  frequency. 
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Abstract — A  linear  theory  using  Laplace  transforms  is  pre¬ 
sented  which  is  applicable  to  both  gyrotron  traveling  wave  am¬ 
plifiers  (gyro-T>VT’s)  and  gyrotron  backward  wave  oscillators 
(gyro-BWO’s).  The  validity  of  the  linear  theory  is  verified  by 
comparing  it  with  an  existing  nonlinear  self-consistent  theory 
based  on  a  different  approach.  In  coidunction  with  a  time- 
dependent  multimode  particle  simulation  code,  the  linear  theory 
is  applied  to  study  the  stability  of  harmonic  gyro-TWT's.  It 
is  shown  that  a  harmonic  gyro-TWT  can  be  made  stable  to 
all  forms  of  spontaneous  oscillations  by  employing  a  multistage 
interaction  structure  and  that  it  can  generate  power  levels  far 
in  excess  of  those  possible  for  a  fundamental  gyro-TWT.  The 
linear  bandwidth  of  a  second-harmonic  gyro-TWT  amplifier  b 
also  calculated. 


I.  Introduction 

Being  an  amplifier  capable  of  high-power  millimeter- 
wave  generation,  the  gyro-TWT  is  finding  an  increasing 
number  of  applications,  particularly  in  radar  and  communi¬ 
cation  systems,  owing  mainly  to  its  promise  of  high  gain 
and  broad  bandwidth.  Yet,  in  contrast  to  the  more  developed 
gyrotron  oscillators,  gyro-TWT  amplifier  research  is  still  in 
a  very  preliminary  state.  To  date,  only  a  few  fundamental 
gyro-TWT  experiments  have  been  completed  [l]-{3],  and  the 
gain,  efficiency,  and  bandwidth  of  these  amplifiers  are  still 
significantly  below  theoretical  predictions.  A  very  important 
obstacle  has  been  the  gyro-TWT’s  susceptibility  to  oscillation. 

TVvo  types  of  oscillations  have  been  observed  to  degrade 
the  performance  experimentally.  In  an  early  TEoi  gyro-TWT 
experiment  at  the  Naval  Research  Laboratory  [4],  oscillation 
near  cutoff  caused  by  the  absolute  instability  [5]  occurred  and 
was  stabilized  in  a  later  experiment  by  using  a  waveguide  with 
a  resistive  wall  [6].  The  second  kind  of  oscillation,  gyrotron 
backward  wave  oscillation  (gyro-BWO),  was  observed  in 
a  recent  gyro-TWT  experiment  performed  at  the  National 
Tsing  Hua  University,  in  Taiwan  [2],  where  clear  evidence 
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of  competition  with  the  amplified  TEn  circular  waveguide 
mode  fi'om  backward  wave  oscillation  in  the  TE21  mode  at  the 
second  harmonic  was  recorded.  The  experiment  also  indicated 
the  effectiveness  of  an  attenuating  sever  in  suppressing  the 
spontaneous  oscillations. 

Many  theoretical  works  |7)-[23)  on  the  general  physics  of 
gyro  devices,  including  gyro-TWT’s,  have  been  published, 
but  only  a  few  theoretical  studies  on  the  problem  of  sponta¬ 
neous  oscillations  have  appeared  [24],  [25).  These  theoretical 
analyses  provide  important  insight  into  the  absolute  oscillation 
threshold  for  gyro-TWT' s.  Meanwhile,  a  linear  theory  has 
been  formed  for  the  design  of  a  gyro-BWO  [26],  where  a 
wave  is  excited  which  propagates  toward  the  electron  gun.  The 
stability  of  a  gyro-TWT  can  also  be  jeopardized  by  gyro-BWO 
interactions  in  many  possible  modes  and  at  various  harmonics. 
The  purpose  of  this  paper  is  to  apply  a  similar  linear  theory  to 
study  the  stability  problem  of  harmonic  gyro-TWT’s  operating 
in  any  TEmn  circular  waveguide  mode.  In  conjunction  with 
a  multimode  particle  simulation  code  [27),  it  will  be  shown 
that  a  harmonic  gyro-TWT  amplifier  can  be  made  stable  to  all 
types  of  spontaneous  oscillations  by  shortening  its  interaction 
length. 

In  a  fundamental  gyro-TWT,  the  amplified  RF  frequency  w 
is  approximately  equal  to  the  electron  cyclotron  frequency, 
ilc  =  eBoIrmoc,  where  Bo  is  the  axial  magnetic  field, 

7  =  U  -  (/?i  +  relativistic  mass 

factor,  P  =  v/c,  V  is  the  electron  velocity,  c  is  the  speed 
of  light  in  vacuum,  and  e  and  m©  arc  the  electron  charge 
and  rest  mass,  respectively.  The  magnetic  field  required  for 
operation  is  directly  proportional  to  the  operating  frequency. 
The  availability  of  magnetic  field  will,  therefore,  practically 
set  an  upper  limit  on  the  frequency  which  can  be  amplified  by 
a  fundamental  gyro-TWT.  However,  harmonic  devices  will 
extend  this  limit,  making  high-frequency  sub-millimeter-wave 
amplification  possible. 

Because  harmonic  interactions  are  in  general  weaker  than 
the  fundamental  interaction,  they  are  more  stable  to  oscil¬ 
lations  and  therefore  allow  much  higher  levels  of  electron 
beam  current  [28].  Thus,  harmonic  gyro-'TWT’s  can  yield, 
in  principle,  significantly  greater  power  than  a  fundamental 
gyro-TWT.  The  electron  beam  geometry  which  yields  the 
highest  gain  for  harmonic  interaction  is  an  axis-encircling 
beam,  but  the  cusp  guns  which  can  create  these  beams  are  still 
undergoing  development  [29]  and  have  not  yet  demonstrated 
the  requisite  parameters.  However,  it  is  shown  in  this  paper 
and  in  the  accompanying  nonlinear  study  [30]  that  a  magnetron 
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injection  gun  (MIG),  which  is  a  mature  technology  (31 J,  is  also 
suitable  for  high-power  harmonic  gyro-TAVT’s.  The  general 
stability  study  presented  in  this  paper  is  appropriate  for  both 
axis-encircling  and  MIG  beam  geometries. 

This  paper  is  organized  as  follows.  In  Section  11,  a  general 
analytical  linear  theory  is  presented.  The  Laplace  transfor¬ 
mation  approach  follows  that  of  Park  el  al.  (26).  While  the 
equations  given  in  (26)  apply  only  to  the  TE„o  rectangular 
and  TEq,.  circular  waveguide  modes,  the  linear  theory  (32| 
described  in  this  paper  can  be  applied  to  any  TEmn  mode 
of  cylindrical  waveguide,  which  is  the  circuit  geometry  most 
commonly  employed  in  gyro  devices.  The  linear  theory  and  the 
time-dependent  multimode  particle  simulation  code  are  then 
used  in  Section  Ill  to  determine  the  conditions  for  stability.  The 
stability  limits  for  the  interaction  length  are  found.  A  reduction 
of  the  interaction  length  can  effectively  stabilize  (28]  ail  forms 
of  spontaneous  oscillations.  Based  on  the  linear  theory,  the 
initial  working  parameter  space  can  be  searched  and  narrowed 
so  that  a  final  set  of  design  parameters  can  be  determined  and 
the  large-signal  performance  of  the  design  can  be  evaluated  by 
self-consistent,  nonlinear  computer  simulations  as  presented  in 
the  accompanying  paper  |30].  Section  IV  concludes  the  paper 
with  a  brief  discussion. 

II.  Analytical  Linear  Theory 

Consider  either  an  annular  or  an  axis-encircling  tenuous 
beam  of  magnetized  electrons  whose  distribution  function  in 
real  and  momentum  space  is  given  by 

/(r,p,t)  =  /o(r,p)  +  /i(r,p,t),  (1) 

with  /o  being  the  unperturbed  (equilibrium)  and  /i  the  per¬ 
turbed  distribution  caused  by  a  TEmn  circular  waveguide 
mode  written  as 

=  ■^F(z)e^  X  Vx<'mn{T-)  , 

E,=  0,  (2) 

=  J^(^)Vx«'mn(r) 


and 

B,  =  -j^F(z)'I'm„(r)  . 


A  time  dependence  of  the  form  has  been  assumed  with 


and 


^'mn(r)  =  Crr,„Jm(kmnr)e  , 

1 


Cmn  — 


F(Z)^  -J 


1  dF 


„  .a  .  1  a 

Vx  —  -1-  Cff-  — 

dr  r  ao 


where  k  =  w/c,  kmn  =  x^nlfw  and  imn  is  the  nth  root 
such  that  Jm(^mn^tn)  =  0  is  satisfied  at  the  waveguide  wall 
radius,  r^.  For  a  tenuous  electron  beam,  the  transverse  field 
profiles  of  the  waveguide  modes  are  assumed  to  be  the  same 


Fig.  l  Projection  of  equilibrium  eteciron  oibit  on  cross-sectional  ptane  of 
waveguide  in  uniform  magnetic  field  Do  Point  O  is  the  center  of  waveguide, 
r„.  is  the  waveguide  radius,  point  Q  is  the  instantaneous  position  of  the 
electron,  rx.  r^,  and  Gr  are  the  electron  Larmor  radius,  gliding  center  radius, 
and  azimuthal  angle,  respectively,  and  r  and  6  are  real-space^lar  coordinates. 

as  those  of  a  cold  waveguide.  The  function,  F(z),  having  the 
dimension  of  field  strength,  is  the  axial  profile  of  the  mode. 
Consistent  with  a  linear  analysis,  the  assumption  of  /i  «:  /o 
is  also  made. 

Using  the  TEmr>  mode  projection  operator  and 

Graff’s  addition  theorem  for  Bessel  functions,  from  Maxwell’s 
wave  equation  one  can  straightforwardly  obtain 

(£  +  (3) 


with  the  source  term 


Si,)  , 


,=  -=o  -'o 

’  ^  ^T£fryx(r’,  8(^Tnnf'c 


(4) 


where  =  y/{iJJcj^—k^  is  the  axial  wave  number  of 
the  wave,  and  are  the  beam’s  guiding  center  and  Larmor 
radii,  respectively,  Jx{r.  t)  is  the  transverse  component  of  the 
beam  current  density,  and  A  =  uj  t  -  m  rj!  +  (s  —  m)(4>  -  rr) 
is  a  phase  variable  defined  with  reference  to  Fig.  1.  Equation 
(3)  is  the  reduced  wave  equation  for  the  axial  field  profile 
F(z)  driven  by  the  induced  beam  current. 

The  induced  beam  current  resulting  from  the  influence  of 
the  RF  field  is  described  by  the  perturbed  distribution  function 
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/i(>'iP'0.  which  can  be  written  as 

/i(r,p,0  =  ^  f  dz'T,{z  ~  z')F,{z') 

(5) 

when  the  linearized  Vlasov  equation  is  solved  by  the  method 
of  characteristics.  In  (5)  the  transfer  function  T,  is  defined  as 

T,{z  —  z')  =  — ~  ygv 

and  the  effects  of  the  RF  fields  are  represented  by  (7),  given 
at  the  bottom  of  the  page.  Equation  (5)  is  a  typical  hysteresis 
integral  and  indicates  that  the  RF  electric  field  F,  at  a  point 
z'  affects  the  electron  distribution  at  a  point  z  cumulatively 
through  the  transfer  function  T,(2  -  z').  As  a  result  of  the 
modulation  produced  by  the  RF  fields,  the  electron  energies 
are  perturbed  and,  therefore,  so  are  their  cyclotron  frequencies, 
which  induces  an  azimuthal  bunching  in  the  electron  phase 
distribution.  This  electron  phase  modulation  represents  an 
induced  beam  current,  which  is  determined  by  the 

electron  distribution  function  /i(r,p,f)  given  by  (5).  When 
•^x(r,  t)  is  substituted  into  (3)  and  (4),  the  field  amplitude  F{z) 
can  be  determined  self-consistently  by  solving  the  resulting 
one-dimensional  integrodifferential  equation. 

The  integrodifferential  equation  can  be  solved  by  using  the 
Laplace  transformation  defined  as 


rp.dp^Tdp. 

moc^  Jo  Jo  J-ac 

fo  —  (x  ^11^  Ff  ni  (t.'mnrc,  i-mn  rZ.  ) 

”  ^11  )  '^3  m(^*pnn^c»  rz.)  -  ^mn  *  *m  (  n  »  km„TL) 

^  -  kr.  -  ^ 

(H) 


So(fc|i)  =j 


Jo  J — oo  ^ 

(UZ  -  fciiu,  - 

ri) 

w  -  fciitz,  - 


C[F{z)]  =  F(i:||)  ^  jT"  d2e^*MF(^)  (g)  r.„(i,y)  =2H.„(i,y)  +  yj;(y)|2j2_„(a:)j;(y) 

I  is  a  complex  wave  number.  Applying  the  Laplace  -  7  fi'l  + 

lation  and  using  the  convolution  theorem,  (3)  can  be  »vyj  - mt  ;  ,  _ mt  ^  - m 


where  fe))  is  a  complex  wave  number.  Applying  the  Laplace 
transformation  and  using  the  convolution  theorem,  (3)  can  be 
reduced  to  a  simple  algebraic  equation  for  F(k\^)  in  fc||  space 

-  '=11  -  fcL)F(fc||)  =  S(fc,|)  -  jfc||F(0)  +  F'(O)  (9) 


F'{0)  = 


dF{z)\ 


are  the  initial  conditions. 

Under  the  Laplace  transformation,  the  source  term  S{z)  on 
the  right-hand  side  of  (3)  becomes 

S(*„)  =  Si(fc„)f’(fc„)  -f  So(fc„)F(0).  (10) 


l^.m(i,y)  =  -  ^y>^i(y){>/.-i(y)l-/J-m-i(a:) 

-  ^?-n.(x)] 

+  ^.  +  i(y)[^L,„  +  i(x)  -  J^,„(x)]}. 

The  evaluation  of  the  remaining  integrals  depends  upon  the 
specific  functional  form  of  the  electron  beam  distribution 
function  /o(r’c,Px,p«).  For  a  monocnergetic  cold  beam,  the 
initial  distribution  function  can  be  written  as 

fo  =  Nk6{'r-’yo):^ — 6{Px-Pi.o)x^Hrc-rco)  (14) 
jTrpxo  ZitVco 


F.(z')  =J„_,{krn„rMk^nrL)\(u,F+-^kk^„Ap^- 

LV  7mo  J  dp±  7mo  9p*J 

1  ^  Jo 

('^m  — »  — l(^mn^c)«^#+l(fcfnn^Zp)  *An— ■+!  (^mn^c  )*^j- 1  (^mn^Z.))  /  — 
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and  ihe  integrations  in  (11)  and  (12)  are  trivially  removed, 
with  the  integration  variables  being  replaced  by  thei-^  initial 
values.  In  (14),  the  constant  Ni,  is  the  number  of  electrons  in 
one  unit  of  axial  length. 

Making  use  of  the  An  space  representation  of  the  field 
amplitude  /’(A||),  as  determined  by  (9),  of  the  inverse  Laplace 
transformation 

A  1  /•jA  +  OO 

F{z)  =  £-i(F(A||)]  =  —  [  dA||e--’'^ii^F(A|l), 

and  of  the  residue  theorem,  the  field  amplitude  in  real  space 
can  be  expressed  in  terms  of  the  initial  values  F(0)  and  F'(0) 
as 


F(2) 

+  F'(0)y  - 


where 


(15) 


^(*11 )  -  2  -  •5l(A||)  , 


(16) 


“''II 

and 


^(^11 )  =  50(^11)  -  jA||  .  (17) 

In  (15),  the  Aj|i  is  the  tth  zero  of  first  order  of  the  equation 


^(^Ib)  =  0,  (18) 


which  is  simply  the  complex  dispersion  relation  describing  the 
coupling  between  the  electron  beam  and  the  waveguide  modes. 

Given  a  set  of  beam  and  wave  parameters,  the  wave 
evolution  is  completely  determined  by  (18)  and  (15).  For  an 
ideal  beam  with  the  distribution  function  given  by  (14),  the 
dispersion  relation  (eq.  (18))  becomes 


D(u;,fc||)=  -j  ~  k. 


-  k 


2 

mn 


^h/U 

- 1 

1 

0 

1  ^*m(^mn^cOi  ^mnT£,o) 

(w  -  A||1J||0  -  sflco/To)^ 

^  (“'  -  fell^'llo)  'FtmjkfnnrcO,  Lnnr'Lo) 
W  -  A:||V||0  -  sClco/'yo 


_  ^mnV±0  f^»m(l>TT»nr'cOi  Amn^Lo) 


W  -  fcllWllo  -  sflco/^o 


=  0 


(19) 


where  =  (3  x  109)-‘ (moc3/e)/3||o7o  %  kA 

is  the  Alfv6n  critical  current,  /(,  is  the  dc  beam  current 
in  amperes  and  the  constant  =  Jl{x„n){l  - 

^^/^mn)-  Equation  (19)  is  the  same  as  the  previous  result 
[25]  derived  from  a  different  approach.  In  the  above,  it 
has  been  assumed  that  the  axial  magnetic  field  is  tuned 


such  that  the  TE^n  wave  interacts  primarily  at  the  sth 
harmonic  cyclotron  frequency.  It  then  follows  that  the 
term  in  (19)  will  dominate  over  terms  being  proportional  to 
F,„,  and  U,m,  as  it  is  inversely  proportional  to  the  second 
order  of  {w  -  A||U||o  -  silco/lo)  ~  ^1-  'he  coupling 
function  H tmik^n^c,  k^nV i)  that  gives  rise  to  the  cyclotron 
instability.  The  dispersion  relation  (eq.  (19))  is  a  fourth-order 
polynomial  in  ut  and  Ay,  which  has  four  complex  roots.  When 
these  four  roots  are  substituted  into  (15),  a  superimposed  wave 
made  up  of  four  individual  waves  is  formed.  For  the  situation 
where  the  beam  cyclotron  mode  couples  with  the  waveguide 
mode  in  the  forward  direction,  there  exist  one  backward  wave 
with  constant  amplitude  and  three  forward  waves,  of  which 
one  is  growing,  one  is  decaying,  and  the  last  one  remains 
constant  in  amplitude.  The  situation  will  be  reversed  if  the 
beam  cyclotron  mode  couples  with  the  waveguide  mode  in 
the  backward  direction. 

Once  the  wave  amplitude  profile  has  been  determined,  the 
time-averaged  power  flow  in  the  circuit  can  be  calculated  as 


Fu,  =  Rel^j  dA  ■  ^{E  X  F*)| 

-^{F{z)3F'[zr}.> 


1 


S’T  wA^„ 

The  linear  gain  is  then  given  by 

P^{z)  Im{F(2)F'(2)-} 


G(2)  = 


F^(0)  Im{F(0)F'(0)-) 


(20) 


(21) 


for  a  forward  growing  wave  (Re{A|i}  >  0),  and 


G{z) 


P^(0)  Im{F(0)F'(0)*} 

P^z)  Im{F(2)F'(2)*} 


(22) 


for  a  backward  growing  wave  (Re{A||}  <  0).  : 

Equations  (15),  (19),  (21),  and/or  (22)  are  solved  subject  to  j 
the  following  boundary  conditions: 

1)  Forward  Traveling  Wave  Launched  at  z  =  0: 

|F(0)|  = 

and 

-j^lz=o  =  -jkzFiO). 

2)  Backward  Traveling  Wfave  Initialized  at  z  =  L: 


and 


dF{z) 
dz  '*=« 


7A,F(0) 


Re{F(2)}  +  jIm{F(2)}  =0  at  2  =  L.  (23) 


It  is  straightforward  to  find  the  linear  gain  for  the  forward  wave  ; 
by  simply  following  the  wave  evolution  through  the  interaction 
tube.  For  the  backward  wave,  however,  one  must  solve  the 
set  of  simultaneous  nonlinear  equations  (eq.  (23))  for  the  two 
unknowns,  namely,  the  backward  wave  oscillation  frequency 
and  the  start  oscillation  length.  When  both  the  real  and  the  ; 
imaginary  part  of  (23)  are  satisfied,  the  gain  given  by  (22) 
will  be  infinite,  which  implies  that  the  device  will  oscillate, 
yielding  finite  output  power  at  c  =  0. 
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III.  OsciL  .ilON  Conditions 

A  study  of  the  dispersion  relation  reveals  that  two  types 
of  spontaneous  oscillations  may  be  excited  in  a  gyro-TWT 
amplifier,  oscillations  arising  from  the  absolute  instability  and 
gyrotron  backward  wave  oscillations  (gyro-BWO).  In  Fig.  2, 
the  cold  waveguide  dispersion  relation  and  the  beam-wave 
resonance  lines  are  shown  for  both  a  TEj^  and  a 
gyro-TWT  amplifier  (the  superscripts  represent  the  cyclotron 
harmonic  number).  A  100  kV  beam  with  a  =  vx/t^H  =  1 
and  Bo  —  0.98Bg  have  been  assumed,  where  Bg  is  the 
grazing  magnetic  field  at  which  the  beam-wave  resonance  line 
grazes  the  operating  waveguide  mode  curve.  The  solid  curves 
and  solid  straight  lines  represent  the  operating  waveguide 
mode  and  the  cyclotron  harmonic  resonance,  respectively. 
Also  shown  in  the  figure  by  the  broken  curves  and  lines  are 
other  neighboring  modes  and  harmonics.  It  is  seen  that  in  a 
TEji  gyro-TWT,  no  oscillation  at  the  fundamental  frequency 
can  be  excited,  since  the  fundamental  beam-wave  resonance 
line  is  well  below  the  lowest  order  waveguide  mode. 

Wave  amplification  will  occur  when  the  beam  cyclotron 
mode  couples  with  the  waveguide  mode  near  the  grazing 
intersection  point.  There  will  be  a  finite  width  of  unstable 
spectrum  about  the  intersection.  If  this  unstable  band  is  con¬ 
fined  to  the  forward  wave  region  only  (/:||  >  0),  a  wave  will 
travel  down  the  interaction  tube  with  convectively  growing 
amplitude.  However,  if  the  unstable  band  is  sufficiently  wide 
to  extend  into  the  backward  wave  region  (/;||  <  0),  waves 
growing  in  the  backward  direction  will  be  excited  and  provide 
internal  feedback  to  yield  self-oscillation  near  cutoff.  If  this 
oscillation  occurs,  it  will  grow  in  time  (absolutely)  when 
viewed  at  a  fixed  spatial  position. 

Relying  on  the  original  criteria  in  [5],  one  may  distinguish 
the  absolute  instability  from  the  convective  instability.  It  is 
known  that  the  absolute  instability  occurs  as  a  zero  of  second 
order  of  the  dispersion  relation  (eq.  (19)).  When  the  dispersion 
relation  is  solved  and  the  behavior  of  its  roots  is  traced  numer¬ 
ically  in  the  u>  —  A;||  space,  the  threshold  current  for  the  onset 
of  the  absolute  instability  can  be  determined  [25].  Fig.  3  shows 
the  threshold  currents  for  the  absolute  instability  in  (a) 
and  (b)  TEji^  gyro-TWT  amplifiers  as  a  function  of  a.  A  100 
kV,  axis-encircling  beam  was  assumed  in  both  calculations. 
For  a  MIG-type  beam,  the  start  oscillation  currents  will  be 
somewhat  higher  than  those  shown  in  Fig.  3,  as  the  interaction 
between  off-axis  electrons  and  the  wave  is  weaker.  It  is  clear 
from  the  figure  that  the  start  oscillation  current  has  a  strong 
functional  dependence  on  a.  As  the  beam  current  must  be 
maintained  below  the  relatively  low  threshold  current  level  to 
ensure  stability,  the  output  power  capability  of  a  fundamental 
gyro-TWT  is  greatly  limited.  However,  a  harmonic  gyro- 
TWT  has  a  weaker  interaction  strength  and  can  therefore  be 
operated  with  a  much  higher  value  of  current,  as  shown  in 
Fig.  3  (more  than  a  fivefold  increase  in  this  case).  In  other 
words,  a  harmonic  gyro-TWT  amplifier  is  more  stable  to  self- 
oscillation  owing  to  the  absolute  instability  than  a  fundamental 
gyro-TWT.  As  a  result,  a  harmonic  gyro-TWT  can  generate 
significantly  higher  output  power  than  a  fundamental  gyro- 
TWT  for  a  given  beam  voltage  and  velocity  ratio  q. 


Fig.  2.  Uncoupled  dispersion  relation  of  operating  mode  (intersection  of 
unbroken  curves)  and  likely  oscillating  modes  (intersections  of  broken  curves 
with  negaUve  ik|| )  for  (a)  TE]','  gyro-TWT  amplifier  and  (b)  gyro-TWT 

amplifier  (100  kV,  a  =  l,5o  =  O.SSF*). 

In  addition  to  the  self-excited  oscillations  caused  by  the 
absolute  instability,  a  gyro-TWT  amplifier  is  also  susceptible 
to  spontaneous  backward  wave  oscillations  which  occur  at 
inteisections  with  negative  fc||.  As  shown  in  Fig.  2(b),  there 
are  various  gyro-BWO  modes,  TE^i^TEj^^  etc.,  competing 
with  the  operating  mode  in  the  second-harmonic  gyro- 

TWT.  These  modes  will  be  excited  with  the  gain  profiles  given 
by  (22)  if  the  interaction  tube  length  is  longer  than  the  critical 
lengths  determined  by  the  gyro-BWO  boundary  conditions  (eq. 
(23)). 

The  starting  oscillation  length  of  each  competing  gyro-BWO 
mode  has  been  calculated  and  the  results  for  the  TE^j'  and 
TE^,^^  modes,  which  have  the  shortest  critical  lengths,  are 
shown  in  Fig.  4  as  a  function  of  beam  current  for  a  beam  of  100 
kV,  0=1,  rc/r^  =  0.4,  r„,  =  0.44  cm,  and  Bo/Bg  =  0.98. 
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Fig  3  Ocpcndcnct  of  sun  oscillition  current  on  velocity  retto  o  for 
ibsolute  insubiliiy  ui  («)  fundemenul  TEi  1  and  (b)  second -harmon.c  TE,, 
gvro-TVNT  amplifiers  for  two  values  of  magnetic  field  (100  kV,  axis-encircline 
beam)  * 

Notice  from  Fig,  2  that  the  TEj^'  intersection  actually  occurs 
for  a  forward  wave.  Since  the  interaction  is  so  close  to 
the  cutoff,  spontaneous  oscillations  arising  from  the  absolute 
instability  will  almost  definitely  occur  unless  proper  meth¬ 
ods  for  its  suppression  are  used.  Using  the  time-dependent, 
multimode  particle  simulation  code,  we  have  observed  that 
there  also  exists  a  critical  length  for  oscillation  owing  to  the 
absolute  instability.  Since  the  present  analytical  linear  theory 
model  does  not  include  finite  length  effects  for  the  spontaneous 
oscillations  caused  by  the  absolute  instability,  the  magnetic 
field  ratio  Bq/B,  was  artificially  increased  17%.  This  was 
done  so  that  the  TEjj^  intersection  occurs  for  a  backward  wave 
such  that  the  preceding  analytical  theory  can  be  used  to  obtain 
®  "^31  critical  length  comparable  to  that  predicted  by  the 
more  accurate  multimode  particle  simulation  code.  The  critical 
length  of  12  cm  found  from  the  multimode  simulation  code 
for  TE3J  oscillation  for  the  correct  design  value  of  magnetic 
field  has  also  been  included  in  Fig.  4  as  a  solid  circle.  Fig. 
4  shows  that  it  is  the  critical  length  for  the  TEj^’  mode  that 
limits  the  interaction  length.  From  Fig.  4  it  can  be  seen  that 
the  critical  length  is  very  sensitive  to  the  beam  current  in  the 
low-current  regime  (7,^  However,  in  the  high-beam- 

curreni  regime  of  interest,  the  critic,!!  length  depends  onlv  very 


1,(A) 

Fig.  4.  Dependence  of  critical  length  on  beam  current  for  backward  wave 
oscillations  at  second  harmonic  in  TE]  1  mode  and  third  harmonic  in  TE31 
mode  (too  kV,  o  =  =  0.4.  r„,  =  0.44  cm,  BojBg  —  0.98  for 

TE^,^*  plot  and  BajBg  =  1. IS  for  plot).  Solid  dot  is  from  multitiKxle 

partide  simulation. 

weakly  on  the  beam  current.  In  this  regime,  the  interaction 
length  needs  be  reduced  only  slightly  to  maintain  stability  for 
a  very  large  increase  in  the  beam  current. 

Using  an  interaction  length  less  than  the  shortest  critical 
length  for  all  possible  spontaneous  oscillations,  a  gyro-TWT 
amplifier  can  kept  stable.  If  the  critical  length  is  too  short 
to  allow  amplification  with  reasonable  gain,  multiple  stages 
separated  by  attenuating  severs  can  be  used,  where  the  length 
of  each  stage  is  still  kept  shorter  than  the  critical  length.  This 
procedure  (28]  will  ensure  single-mode  operation  together  with 
a  high  amplification  gain. 

IV.  DiscxissioN  AND  Conclusions 

The  prediction  of  the  linear  theory  described  in  Section  II 
has  been  compared  with  that  of  an  existing  self-consistent 
nonlinear  theory  (30).  The  results  arc  shown  in  Fig.  5,  where 
the  power  evolution  in  a  TEj^^  gyro-TWT  amplifier  with 
Bo  =  0.98B,  is  depicted  for  a  100  kV,  25  A  beam  with 
Q  =  1  and  Tc/tu,  =  0.4.  The  comparison  made  in  Fig.  5  shows 
perfect  agreement  in  the  linear  region  between  the  two  theories 
derived  from  completely  different  approaches,  an  indication  of 
the  validity  of  both  theories  in  that  region.  As  expected,  the 
linear  theory  fails  in  the  saturation  region,  where  nonlinearities 
become  important. 

In  addition  to  the  stability  considerations,  the  linear  theory 
has  also  been  used  to  evaluate  the  linear  performance  of  a 
gyro-TWT  amplifier  for  various  beam  parameters.  The 
linear  growth  rate  is  shown  in  Fig.  6  with  axial  magnetic 
field,  beam  current,  and  clcaron  velocity  ratio  as  varying 
parameters.  As  expected,  the  interaction  strength  becomes 
stronger  as  q  or  the  beam  current  is  increased.  The  linear 
growth  rate  also  increases  with  magnetic  field  strength  as  Bo 
is  raised  toward  the  grazing  magnetic  field,  B,.  For  a  25  A,  100 
kV  beam  with  o  =1  and  r^/r^.  =  0.4  interacting  with  a  TE^^* 
mode  wave  in  a  uniform  magnetic  field  of  Bo  ~  0.98Bj,  the 
3  dB  bandwidth  predicted  by  the  linear  theory  is  10.6%, 
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Kig.  5  Predicted  spatial  power  profile  of  35  GHz  wave  in  gyro-TWT 

Irom  linear  theory  (broken  curve)  and  nonlinear  simulation  (solid  curve)  (100 
kV,  25  A,  Q  =  l.rc/r^,  =  0.4.  and  Bo  =  0.98B,). 


Fig  6.  Linear  growth  rate  as  a  function  of  frequency  in  100  kV,  TEI?* 
gyro-TWT  amplifier  for  (a)  /(,  =  25  A  and  o  =  1  for  several  values  of 
magnetic  field,  (b)  /t  =  25  A  and  Bo  =  0.98Bj  for  several  values  of  o. 
and  (c)  o  =  1  and  Bo  =  0.98Bj  for  several  values  of  beam  current.  Ac  and 
■  V  are  the  cutoff  wavelength  and  angular  frequency,  respectively.  The  ratio 
"f  the  beam  guiding  center  radius  to  the  waveguide  wall  radius  is  0.4. 


The  analytical  linear  theory  presented  in  Section  II  has  been 
developed  based  on  a  single-mode  as.sumption  and  employed 
to  analyze  gyro-TWT  stability,  which  is  inherently  a  multi¬ 


Ihl 

mode  problem.  However,  such  a  single-mode  assumption  and 
analysis  are  still  proper  for  stability  considerations,  as  the 
critical  lengths  so  determined  are  less  than  the  start  oscillation 
lengths  for  the  driven  system.  Using  these  critical  lengths, 
each  subsection  of  the  gyro-TWT  will,  ihcrefore,  be  stable  to 
spontaneous  oscillation,  and  the  single-mode  analysis  for  each 
stage  will  then  be  valid. 

A  general  analytical  linear  theory  has  been  developed  for 
gyro  devices  involving  any  TE„,„  circular  waveguide  mode 
at  any  cyclotron  harmonic  number.  It  has  been  demonstrated 
generally  that  harmonic  gyro-TWT  amplifiers  not  only  can 
amplify  high  frequency  at  reduced  magnetic  field  with  ample 
gain,  but  also  have  potential  for  generating  extremely  high 
power  in  the  millimeter-wave  region.  A  second-harmonic 
gyro-TWT  amplifier  can  be  kept  stable  for  beam  current  levels 
one  order  of  magnitude  higher  than  levels  permitted  in  a  stable 
fundamental  gyro-TWT.  The  analytical  linear  theory  can  give 
a  thorough  physical  understanding  of  the  mechanisms  of  the 
various  spontaneous  oscillations  and  forms  the  basis  for  the 
design  of  high-power  harmonic  gyro-TWT  amplifiers. 

• 
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High-Power  Harmonic  Gyro-TWT’s — Part  II: 
Nonlinear  Theory  and  Design 
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Abstract —  Based  upon  an  analytical  study  of  the  stability 
problems  of  gyrotron  traveling  wave  amplifiers  (gy^o-TV^T’s), 
the  design  of  an  extremely  high  power  second-harmonic  gyro- 
TWT  has  been  developed  and  then  evaluated  and  optimized  with 
a  self-consistent  nonlinear  numerical  simulation  code.  Although 
an  axis-encircling  electron  beam  yields  stronger  interaction,  a 
design  based  on  the  more  common  magnetron  injection  gun 
(MIG)  type  beam  is  presented.  Using  a  100  kV,  25  A  MIG  beam 
with  a  =  1  and  an  axial  velocity  spread  of  S%,  nonlinear  self- 
consistent  analysis  of  a  three-stage  second-harmonic  gyro-TWT 
amplifier  predicts  a  peak  output  power  of 533  kW,  peak  efficiency 
of  2U%  and  a  7.4%  saturated  bandwidth,  which  verifies  the 
assessments  of  analyticai  theory  that  a  stable  harmonic  gyro- 
TWT  can  generate  power  levels  an  order  of  magnitude  higher 
than  those  possible  from  a  fundamental  gyro-TWT.  It  is  shown 
that  the  positioning  of  the  eiectron  beam  is  very  important.  A 
multistage  structure  is  employed  in  the  design  to  recover  the  loss 
in  gain  resulting  from  the  shortening  of  the  interaction  sections 
to  ensure  stability. 

I.  Introduction 

The  gyrotron  traveling  wave  amplifier  (gyro-TWT)  has 
been  shown  both  theoretically  and  experimentally  to  be  a 
highly  efficient  millimeter-wave  amplifier.  Because  it  employs 
a  simple  smooth  waveguide  structure,  the  gyro-TWT  has  far 
greater  power  handling  capability  than  conventional  TWT’s. 
Having  substantial  gain  and  a  broad  bandwidth,  gyro-TWT’s 
have  found  a  wide  variety  of  applications  such  as  power  ampli¬ 
fiers  in  millimeter-wave  radar  and  communication  systems.  To 
dale,  several  experimental  gyro-TWT’s  have  been  reported.  An 
early  gyro-TWT  developed  at  the  Naval  Research  Laboratory 
(NRL)  (l]-[4]  has  amplified  a  TEqi  circular  waveguide  mode 
near  35  GHz  with  17-32  dB  gain  over  a  13%  bandwidth  with  a 
maximum  efficiency  of  16%.  A  recent  gyro-TWT  experiment 
at  NRL  f5j  demonstrated  the  capability  of  broad  bandwidth 
(33%  bandwidth  with  20  dB  linear  gain  at  32.5  GHz)  by 
tapering  both  the  axial  magnetic  field  and  the  waveguide.  The 
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first  Varian  gyro-TWT  [61-[81  achieved  a  remarkable  120  kW' 
peak  power  near  5  GHz  using  the  TEn  circular  waveguide 
mode  with  a  maximum  efficiency  of  26%  and  a  bandwidth 
of  6%.  The  second  Varian  gyro-TWT  [9]  generated  30  kW  at 
94  GHz  with  a  gain  of  30  dB.  Most  recently,  a  group  at  the 
National  Tsing  Hua  University  has  investigated  the  stability 
problem  by  demonstrating  a  27  kW,  K,-band,  severed  gyro- 
TWT  amplifier  with  16%  efficiency,  35  dB  gain,  and  a  full 
width  at  half  maximum  (FWHM)  bandwidth  of  7.5%  using 
the  TEii  circular  waveguide  mode  [10].  All  these  reported 
gyro-TWT’s  have  produced  moderately  high  levels  of  output 
power  (10-100  kW)  and  involved  amplification  only  at  the 
fundamental  cyclotron  harmonic  frequency.  ' 

High-harmonic  devices  have  also  been  developed  to  allow 
high-freque ;icy  operation  at  reduced  magnetic  field  [llj-{16]. 
These  high-harmonic  devices  used  large-orbit,  axis-encircling 
electron  beams  for  strong  interaction.  However,  the  technique 
for  forming  such  beams  at  high  levels  of  current  with  the  small 
axial  velocity  spreads  required  by  an  amplifier  is  still  under 
development. 

Important  for  a  gyro-TWT’s  stability,  the  working  beam 
current  must  be  limited  below  the  threshold  value  for  an  ab¬ 
solute  instability  to  occur;  also,  the  interaction  length  must  be 
kept  shorter  than  a  critical  length  beyond  which  spontaneous 
gyrotron  backward  wave  oscillations  will  take  place  [17].  In 
an  sth-harmonic  cylindrical  TE^i  gyro-TWT,  first-harmonic 
interaction  will  not  occur,  since  the  fundamental  cyclotron 
resonance  line  is  well  below  the  lowest  order  waveguide 
mode  as  viewed  in  the  uncoupled  dispersion  diagram.  Previous 
theoretical  studies  have  concluded  that  a  harmonic  gyro-TWT 
is  more  stable  to  oscillation  than  a  fundamental  gyro-TWT 
amplifier  [17],  [18],  since  harmonic  interactions  are  in  general 
weaker  than  the  fundamental  interaction.  Therefore,  harmonic 
interactions  allow  significantly  higher  levels  of  beam  current 
and  output  power.  If  the  gain  of  a  single-stage  gyro-TWT 
with  a  limited  interaction  length  is  too  low,  multiple  stages 
separated  by  attenuating  severs  can  be  used  to  boost  the 
amplification  gain.  Based  on  this  marginal  stability  design 
procedure,  a  high-power,  second- harmonic  gyro-TWT  design 
using  an  axis-encircling  beam  is  given  in  [17].  Although 
an  axis-encircling  beam  geometry  is  preferable  for  stronger 
interaction,  cusp  guns  which  can  create  these  beams  are  still  in 
development  [19]  and  have  not  yet  demonstrated  the  requisite 
parameters.  Fortunately,  nearly  equivalent  performance  levels 
at  the  second  harmonic  can  be  achieved  utilizing  an  annular 
beam  formed  bv  a  conventional  MIG.  Since  MiG's  are  a 
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mature  technology  [2()j,  the  high-power,  second-harmonic 
gyro-TWT  proposed  in  this  paper  employs  a  MIG  as  its 
electron  gun. 

This  paper  is  organized  as  follows.  In  Section  II,  a  presenta¬ 
tion  of  the  general  three-dimensional,  nonlinear  self-consistent 
theory  will  be  given  only  briefly,  since  many  nonlinear  theories 
have  appeared  in  the  literature  (21j-{24],  The  slow  time-scale 
approach  followed  in  developing  the  nonlinear  formulation  is 
similar  to  those  used  in  (21)  and  (22).  As  in  (23)  and  (22), 
the  formulation  includes  guiding  center  motion  and  spread, 
an  axial  magnetic  field  taper,  and  the  modeling  of  attenuating 
severs  in  a  smooth-wall  circular  waveguide.  The  design  of 
a  second-harmonic  gyro-TWT  will  be  described  in  Section 
Ill.  Using  the  large  signal  code,  it  is  shown  that  an  output 
power  far  greater  than  that  of  state-of-the-art  fundamental 
gyro-TWT’s  can  be  achieved  for  the  designed  three-stage 
gyro-TWT  amplifier  utilizing  a  MIG  beam.  Finally,  Section  IV 
will  conclude  the  paper  with  a  brief  discussion  of  the  results 
and  possible  future  efforts. 

II.  Nonunear  Theory 

Consider  an  annular  or  solid  electron  beam  interacting  with 
a  copropagating  circularly  polarized  circular  waveguide  TEmn 
mode  wave.  The  electron  beam  is  injected  into  the  interaction 
tube  under  the  guidance  of  an  applied  dc  magnetic  field 
Bq.  The  electrons,  moving  along  helical  trajectories  about 
guiding  centers  located  at  a  radius  r^o,  have  a  substantial 
pan  of  their  kinetic  energy  in  the  form  of  transverse  cyclotron 
motion  and  the  rest  in  the  form  of  axial  motion.  Our  nonlinear 
model  assumes  that  all  electrons  have  the  same  kinetic  energy, 
moc^(7  -  1),  but  the  panition  between  the  transverse 

and  longitudinal  velocities  may  be  different  from  elearon  to 
electron,  where  mo  is  the  electron  rest  mass,  7  =  [1  -  + 

is  the  relativistic  mass  factor,  0  =  v/c,v  is  the 
electron  velocity,  and  c  is  the  speed  of  light  in  vacuum.  It 
is  also  assumed  that  the  electron  beam  is  sufficiently  tenuous 
that  it  does  not  modify  the  normal  mode  of  the  cold  waveguide 
and  its  space-charge  field  can  be  neglected.  Only  a  single  RF 
mode  is  assumed  to  be  present  in  the  waveguide  with  the 
specified  frequency  w.  Attenuation  in  the  waveguide  wall  has 
been  neglected. 

For  the  purpose  of  frequency  scaling  and  also  to  simplify  the 
appearance  of  the  set  of  nonlinear  equations  to  be  presented 
below,  a  normalization  scheme  is  employed  where  the  lengths 
are  scaled  against  some  arbitrarily  specified  length  b  (which 
is  commonly  chosen  to  be  the  waveguide  radius,  r^).  All 
normalized  symbols  are  denoted  by  a  bar.  The  normalization 
scheme  is  summarized  in  Table  1.  All  dimensionless  quantities 
such  as  o(=  vx/v\\),P,  and  7  are  unchanged. 

The  electron  dynamics  are  governed  by  the  Lorentz  force 
equation,  given  by 


IABLE  I 

Nobmaleation  for  Nonlinear  Formulation 


dt 


-p  X  (B  -F  Bo 


where  the  applied  magnetic  field,  Bq,  can  be  expressed  as 

Bo  =  Bo(i)c,  -  -  — — rcr. 

2  uz 


Lengih 

Fiequenc^ 

Time 

Wave  number 
Momemum 


Cyclotron  frequency 


Beam  current 
Cuneni  density 


r  =  r/i, 
i  =  trfh 

^  rn  !•  “■ 

p  =  p/rjioc  =  iP 
E  =  Eeb/rnoc^ 

(irO  -  ^0i  = 
cB{3tb/mQC^ 

h  =  he/mor^ 

J  =  Jett^ /mor^ 


In  the  above,  E  and  B  are  the  RF  fields  of  the  TEmn  circular 
waveguide  mode  as  specified  by 

E±  =  7 - F{z)Cz  X  Vx’^mn)*'), 

^mn 

Ez=0, 

1  dF  ,  . 


B,  =  -j%^F(z)«l'mn(r).  (3) 

K 

A  time  dependence  of  the  form  has  been  assumed  with 

_ 1 _ 

\/7r(i^„  -  rn^Jm(x„n) 


^  .  d  -  1  5 

Vi  = 

dr  r  do 

where  k  =  w/c,  km„  =  Zmn/r„„  and  imn  is  the  nth  root 
such  that  Jm(kmnru)  =  0  is  satisfied  at  the  waveguide  wall 
radius,  Tu,.  The  function  F(z),  having  the  dimension  of  field 
strength,  is  the  unknown  axial  profile  of  the  TEmn  mode. 

Transformation  of  the  coordinates  from  the  cylindrical  sys¬ 
tem  to  the  guiding  center  system  depicted  in  Fig.  1  and 
application  of  Graff’s  addition  theorem  for  Bessel  functions 
will  change  the  general  Lorentz  force  equation  into  a  set  of 
slow  time-scale  equations  of  motion  governing  the  electron 
momentum  and  the  phase  angle,  which  are  given  by 

1  dfleO  -  j  (i-  f  ) 
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and 
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where  p,  and  p^  are  electron  axial  and  transverse  momenta, 
respectively,  A  —  u!l  —  inx'  +  [s  -  rn)((;i  —  t  )  is  a  phase  variable 
defined  with  reference  to  Fig.  1,  is  the  electron  Larmor 
radius,  and  =  f2ru/"i  =  e/io/^rripf  is  the  relativistic 
electron  cyclotron  frequency.  It  has  been  assumed  that  the 
contribution  from  the  sth  harmonic  to  the  interaction  will 
dominate  all  the  other  harmonics.  Similarly,  the  equations 
for  the  electron  guiding  center  motion  can  be  written  in 
normalized  form  as 


dz 
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In  addition  to  the  Lorentz  force  equation,  it  is  also  useful 
to  consider  the  energy  change  equation 

moc^  ~  =  —ev  ■  E, 
at 

which,  after  similar  manipulation,  can  be  rewritten  as 

0x  -  -  _  ,  - 

^  )  cos  0.  (9) 

Equations  (4)-(9)  describe  the  electron  motion  in  the  presence 
of  a  strong  axial  magnetic  field  and  ihe  RF  fields.  A  self- 
consistent  formulation  is  completed  when  the  amplitude  of  the 
wave  field  is  calculated  from  the  induced  ac  current  density. 
This  is  done  by  lengthy  but  straightforward  manipulations, 
very  similar  to  those  leading  to  (4)  through  (9),  to  rewrite  the 
wave  equation  for  the  axial  component  of  the  magnetic  field  in 
terms  of  the  slow  time-scale  variables.  The  normalized  wave 
equations  can  be  obtained  as  follows; 

id^  ~  kmn  -  k^']f  = 


fUo  f’x 
(f'zo)  ll'.l 


~  si  ^ s  (  k] 


Fig.  1  Projection  of  equilibrium  elcciron  orbit  on  cross-sectional  plane  of 
waveguide  in  uniform  magnetic  field  Bq.  Pome  O  is  center  of  Ihe  waveguide. 
r„.  is  Ihe  waveguide  radius,  and  point  Q  is  the  instantaneous  position  of  Ihe 
electron  r^,  tr.  and  Be  ate  'he  electron  Larmor  radius.  Ihe  guiding  center 
radius,  and  the  azimuthal  angle,  respeclivcly,  and  r  and  fr  are  real-space  polar 
coordinates. 


and 

^(kP)  =87r/|,fc^„wCmn/ 

•  (■r:^—Jm-,{kmr,fc)J',ik„„fL)cos§\ 
\(t'a0)|Va|  / 

(11) 

In  deriving  (10)  and  (11),  which  determine  the  wave’s  am¬ 
plitude  and  phase,  as  well  as  (4)-(9),  which  describe  electron 
motion,  the  normalized  axial  spatial  profile  of  the  RF  wave 
field  has  been  assumed  to  be  of  the  form 

F{z)  =  (12) 


where  (^(z)  =  k^z-ipiz),  k{z)  =  dC,ldz  =  —dififdz,  and  ip 
is  the  wave  phase.  Also  9  =  A-((z)  is  the  electron  phase  with 
respect  to  the  wave  and  (  )  stands  for  the  ensemble  average 
over  the  initial  electron  distributions  of  the  axial  momenta, 
guiding  center  positions,  and  gyro  phases. 

Finally,  the  spatial  rate  of  change  of  power  flow  along 
the  waveguide  is  calculated  from  the  Poynting  theorem  and 
appears  as 

=  !bCmnf(~^^J„^-,(k,„r,r^)J',(k„nfi)COs9 

dz  \(i^«o)k2l 

(13) 

where  (vio)  is  the  electrons’  average  initial  velocity,  and  h 
is  the  dc  beam  current. 


Equations  (4)-(13)  form  a  closed  set  of  ordinary  differential 
equations.  Two  of  them  ((9)  and  (13))  are  auxiliary  equations 
which  are  included  to  monitor  the  accuracy  of  the  numerical 
calculations.  The  nonlinear  behavior  of  the  electron  motion  as 
well  as  the  wave  evolution  in  gyro-T\VT’s,  CARM  amplifiers, 
etc.  is  completely  and  self-consistently  determined  by  solving 
this  set  of  equations  Wc  employ  a  fourth-order  Runpe-Kutta 
algorithm  to  solve  the  differential  equations,  subject  to  the 
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initial  conditions  that  k  —  k.  and  df  /dz  =  0,  as  there  is  no  net 
energy  transfer  for  the  initially  unbunched  beam  assumed  at 
the  input  (e  =  0).  The  wave  amplitude  f{z  =  0)  is  calculated 
from  the  input  signal  power.  The  initial  pha  of  the  wave  can 
be  chosen  arbitrarily,  because  all  relevant  quantities  have  been 
averaged  over  the  initial  electron  phases. 

The  initial  state  of  the  electrons  is  chosen  in  such  a  way  that 
the  injection  of  a  monoenergetic,  axisymmetric  beam  is  well 
represented.  It  is  assumed  that  a  finite  number  of  electrons  are 
uniformly  distributed  in  phase,  Ao,  over  a  Larmor  orbit  with  a 
radius  rio.  They  are  also  uniformly  distributed  over  the  ring 
of  an  average  guiding  center  radius.  There  may  be  a  finite 
amount  of  spread  in  the  guiding  center  radius.  For  an  injected 
warm  beam,  a  Gaussian  distribution  function  of  the  form 

Go(p.)  =  >16(7 

is  assumed  for  the  beam  axial  momentum  distribution.  In  the 
above,  i4  is  a  normalization  constant  and  p,o  is  the  average 
initial  axial  momentum. 

III.  Design  of  a  Stable  0.5  MW 
Second-Harmonic  Gyro-TWT 

Application  of  the  above  nonlinear  theory  to  a  100  kV,  25  A 
beam  with  q  =  1  and  Tc/tu,  =  0.4  in  a  single-stage  second- 
harmonic  TE21  gyro-TWT  amplifier  with  waveguide  radius 
r^.  ~  0.44  cm  and  Bq  =  0.98Bg  is  shown  in  Fig.  2,  where 
the  spatial  power  profile  is  plotted.  Bg  is  the  magnetic  field 
which  results  in  a  grazing  intersection  between  the  operating 
mode  and  the  operating  cyclotron  harmonic  resonance  line. 
Also  included  in  Fig.  2  is  the  result  calculated  from  the  linear 
theory  [18]  for  the  same  set  of  beam  and  wave  parameters. 
It  can  be  seen  that  the  two  theories  formed  from  different 
approaches  agree  very  well  in  the  linear  region,  an  indication 
of  the  validity  of  both  theories  in  that  region.  As  expected, 
when  nonlinearities  become  important  in  the  saturation  region, 
where  the  linear  theory  fails,  an  accurate  nonlinear  theory  is 
needed  to  assess  the  large-signal  behavior. 

As  pointed  out  in  the  introduction,  harmonic  interaction  is 
relatively  weak,  leading  to  lower  gain  and  a  longer  interaction 
length  than  the  conventional  fundamental  gyro-TWT  amplifier, 
unless  other  methods  are  applied.  For  example,  for  the  case 
shown  in  Fig.  2,  where  the  critical  length  for  spontaneous 
oscillations  to  occur  has  been  found  to  be  12  cm  [18],  an  input 
power  of  more  than  30  kW  is  required  to  drive  the  device  into 
saturation,  resulting  in  a  power  gain  of  only  14  dB.  It  is  then 
essential  to  use  a  multistage  gyro-TWT  amplifier  if  high  gain 
is  needed.  An  accurate  simulation  of  this  system  requires  a 
nonlinear  theory  which  is  capable  of  modeling  the  attenuating 
severs  which  divide  the  gyro-TWT  into  multiple  sections. 

To  model  a  severed  gyro-TWT,  the  nonlinear  formulation 
presented  in  Section  II  can  still  be  applied  to  each  individual 
stage.  However,  inside  the  sever  region,  the  formulation  has 
been  modified  to  accommodate  the  assumption  of  infinite 
attenuation  such  that 


Fig.  2.  Predicted  spatial  power  profile  of  35  GHz  wave  in  TE^  I  gyro-TWT 
from  linear  theory  (broken  curve)  and  nonlinear  simulation  (solid  curve)  (or 
input  powers  of  100  W  and  30  kW  (100  kV.  25  A.  o  =  1,  r,/r„.  =  0.4. 
r„.  =  0.44  cm,  and  Bo  =  0  OSB,) 

inside  the  sever.  The  electron  beam  dynamics  are  still  deter¬ 
mined  by  (4)-(9)  except  that  all  the  field  quantities  are  set 
to  zero.  The  assumption  of  infinite  attenuation  is  reasonable 
inside  a  sever,  because  the  presence  of  either  lossy  materials 
or  the  perturbation  to  the  waveguide  geometry  will  destroy 
the  electron-wave  synchronism.  As  a  result,  the  effective  field 
seen  by  the  electrons  will  be  essentially  zero  insider  a  sever 
even  if  the  wave  attenuation  is  finite.  Because  of  the  interaction 
with  the  RF  field  prior  to  the  sever  section,  the  electrons  arc 
azimuthally  bunched  in  phase  and  will  inertially  preserve  the 
bunch  through  the  sever  region. 

At  the  exit  of  the  sever,  the  ac  electron  beam  parameters 
are  continuous.  This  alternating  current  will  then  coherently 
induce  an  RF  field  in  the  circuit  of  the  next  stage.  At  the 
entrance  to  the  next  stage,  it  is  assumed  that  the  amplitude 
of  the  induced  RF  field  is  small  (approaching  zero)  and 
dF/dz  =  -jk^F.  In  the  actual  numerical  simulations,  the 
wave  amplitude  /  has  been  set  to  be  a  very  small  value  (<  —40 
dB  of  the  initial  value  at  z  =  0),  but  still  sufficiently  large  to 
avoid  possible  computer  overflow  difficulty.  The  phase  of  the 
wave  at  the  entrance  is  set  to  the  value  at  which  the  wave 
grows  fastest,  as  the  electron  beam  has  been  phase  bunched 
by  the  RF  field  in  the  previous  stages. 

Using  both  the  linear  theory  developed  in  [18]  and  the 
nonlinear  theory  described  in  the  previous  section,  extensive 
parameter  space  searches  have  been  performed,  leading  to 
the  set  of  design  parameters  listed  in  Table  11,  upon  which 
an  experiment  under  construction  at  UCLA  is  based.  All 
numerical  simulations  to  be  shown  in  this  section  are  for  the 
case  listed  in  the  table  unless  specified  otherwise. 

The  MIG  beam  used  in  the  design  has  a  guiding  center 
radius  =  0.4r„„  which  is  chosen  by  careful  study  of 
the  behavior  of  the  coupling  function  //.,m ( Lmn Tr •  kmnTi) 
defined  as  [25],  [18] 

which  has  the  physical  meaning  ('f  bcini;  the  cfieclivt 
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TABLE  II 

Second-Harmonic  35  GHZ  Gyro-TUT  AMPunER  Design 


Beam  voltage 

100  kV 

Beam  current 

35  A 

Magnetic  field 

b,4  kC. 

Center  frequency 

-35  GHe 

Mode 

TE;, 

Cyclotron  harmonic 

2nd 

Circuit  radius 

0.44  cm 

rc/r„. 

04 

0 

1.0 

B/B, 

0.98 

At.||/u|| 

59f 

Number  of  stages 

3 

Section  length 

1 1.5  cm 

Sever  length 

3  cm 

Circuit  length 

40.5  cm 

Simulation  Results: 

Output  power 

533  kW 

Saturated  gain 

54  dB 

Efficiency 

21.3% 

Constant  drive  bandwidth 

6% 

Saturated  bandwidth 

7.4% 

sth  harmonic  component  of  the  TEmn  mode  electric  field 
seen  by  the  electrons.  Fig.  3  shows  the  contour  plots  of 
fi,m{kmnrc,  k,„„ri)  for  the  operating  mode  and  its 

twO(^strongest  competing  oscillation  modes,  the  and 

TE31 ,  where  the  superscripts  represent  the  cyclotron  harmonic 
numbers.  The  dashed  lines  on  the  figures  represent  the  electron 
Larmor  radius  corresponding  to  an  axial  magnetic  field  Bq  = 
0.98Bj  for  the  operating  TEj*'  mode.  Clearly,  the  positioning 
of  the  electron  guiding  center  is  extremely  important,  since  it 
changes  the  relative  interaction  strength  of  all  the  modes.  A 
properly  placed  beam  will  maximize  the  value  of  for 
the  operating  mode  {H2‘2  in  this  case),  while  minimizing  the 
s  for  all  the  other  modes.  Given  an  electron  Larmor 
radius.  Fig.  3(a)  shows  that  H22  is  maximum  for  an  axis- 
encircling  beam.  For  such  a  beam,  H21  vanishes  while  H33 
is  about  six  times  smaller  than  ^22-  A  TE^^j^  gyro-TWT 
amplifier  design  based  upon  this  beam  geometry  has  recently 
been  developed  with  appreciable  power,  efficiency,  and  gain 
[17J.  However,  difficulties  encountered  in  forming  such  a  high- 
current  axis-encircling  beam  with  low  velocity  spread  still 
need  to  be  overcome. 

As  the  beam  moves  away  from  the  waveguide  center  so  that 
it  corresponds  to  the  more  common  annular  beam  formed  by 
a  MIG,  a  much  more  mature  technology,  both  H22  and  H33 
will  decrease  while  H21  increases.  Notice  that  i/33  decreases 
more  rapidly  than  H22  does,  resulting  in  an  increase  in  the 
ratio  H22/H33.  For  the  chosen  ratio  of  guiding  center  radius 
to  waveguide  wall  radius,  =  0.4,  both  i/33  and  i/21  are 

at  least  one  order  of  magnitude  smaller  than  i/22. 

The  power  evolution  in  the  three-stage  TE^^'  gyro-TWT 
obtained  from  simulation  is  shown  in  Fig.  4  for  three  different 
input  powers  (0.5  W,  1.9  W,  and  4  W).  The  length  of  each 


Fig.  3.  Contour  plot  for  coupling  functions  (a)  //22.  (b)  W33.  »nd  (c) 
Wii.  Absolute  values  of  have  been  magnified  by  100.  Dashed  line 

corresponds  to  100  kV,  o  =  1.  and  Bo/B,  =  0.98. 


section  has  been  chosen  to  be  less  than  the  critical  length 
for  TE3^^  oscillation,  which  has  been  found  to  be  the  most 
threatening  spontaneous  oscillation  [18].  As  in  a  single-stage 
gyro-TWT  amplifier,  the  output  power  saturates  at  shorter 
distance  as  the  input  power  is  increased.  It  is  found  that  the 
electrons  will  suffer  from  overbunching  if  the  input  power 
is  too  strong,  resulting  in  lower  saturated  output  power  [27]. 
Simulation  results  have  indicated  that  it  is  easier  to  overdrive 
the  electron  beam  in  a  multistage  device  because  saturation 
occurs  close  to  the  sever  in  the  last  section.  For  the  case  shown 
in  Fig.  4,  the  electron  beam  is  highly  overdriven  when  the  input 
power  is  increased  merely  3  dB  from  1.9  W  to  4  W. 

Shown  in  Fig.  5  is  the  constant  drive  bandwidth  of  the 
second-harmonic  gyro-TWT  amplifier  for  several  values  of 
axial  velocity  spread  (Avjv^  =  0,  3%,  and  5%),  where  all 
parameters  except  the  RF  input  signal  frequency  arc  held 
constant.  The  input  power  is  set  to  the  value  at  which  the 
amplifier  saturates  at  35  GHz.  Clearly,  the  beam  axial  velocity 
spread  narrows  the  constant  drive  bandwidth.  Whereas  an  ideal 
beam  without  a  velocity  spread  gives  a  FWHM  bandwidth 
of  7.1%,  a  beam  with  an  axial  velocity  spread  of  5%  yields 
a  FWHM  bandwidth  of  5.2%.  Fig.  6  shows  the  predicted 
saturated  bandwidth  of  the  harmonic  gyro-TWT  amplifier 
design  for  several  values  of  axial  velocity  spread.  The  saturated 
bandwidths  were  obtained  for  a  fixed  device  length  (40.5 
cm)  with  varying  input  power  such  that  the  maximum  power 
appeared  at  the  output  end  of  the  amplifier.  Again  the  FWHM 
bandwidth  is  reduced  from  13%  to  7.4%  as  the  axial  velocity 
spread  increases  from  0%  to  5‘^( .  The  velocity  spread  has  its 
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Fig.  4.  Predicted  spatial  power  profile  of  35  GHa  wave  in  proposed 
three-stage  gyro-TWT  (Table  II)  for  input  powers  of  0.5  W  (dashed 

curve),  1.9  W  (solid  curve),  and  4  W  (dotted  curve). 
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Frequency  (  GHz ) 

Fig.  5.  Constant  drive  bandwidth  of  proposed  three-stage  gyro-TWT 

(Table  II)  for  axial  velocity  spreads  of  0%  (dashed  curve),  3%  (dotted  curve), 
and  5%  (solid  curve). 

greatest  effect  in  the  high-frequency  region  of  the  bandwidth. 
This  is  because  the  larger  wave  numbers  associated  with  the 
higher  frequencies  yield  a  larger  Doppler  shift  spread  and 
therefore  a  greater  mismatch  of  the  beam-wave  resonance 
condition.  For  the  case  of  Av^/vz  =  5%,  a  peak  output  power 
of  533  kW  is  predicted  at  34.5  GHz,  which  is  almost  an  order 
of  magnitude  higher  than  the  power  levels  achieved  by  state- 
of-the-art  fundamental  gyro-TWT’s.  At  34.5  GHz,  the  peak 
efficiency  and  gain  arc  21.3%  and  54  dB,  respectively. 

IV.  Discussion  and  Conclusions 

Although  reducing  the  interaction  length  is  a  general,  ef¬ 
fective,  and  practical  way  to  suppress  all  spontaneous  oscilla¬ 
tions  in  a  high-power  gyro-TWT,  external  feedback  caused 
by  reflections  from  poor  input/output  wave  couplers,  poor 
sever  matching  circuits,  etc.,  could  induce  oscillation,  render¬ 
ing  wave  amplification  impossible.  These  external  reflections, 
however,  have  been  neglected  in  the  nonlinear  simulations. 
The  critical  interaction  length  for  oscillation  can  be  expected 
to  decrease  for  nonzero  reflection  coefficient,  but  the  technique 
of  using  multiple  stable  sections  to  boost  the  gain  is  still 


Frequency  ( GHz ) 


Fig  6.  Saturated  bandwidth  of  proposed  three-stage  TE^j*  gyro-TWT 
(Table  II)  for  axial  velocity  spreads  of  0%  (dashed  curve).  3%  (dotted  curve), 
and  5%  (solid  curve). 

appropriate.  To  prevent  oscillations  from  occurring,  the  sum  of 
the  return  loss  from  each  sever  and  its  complement,  a  coupler 
or  another  sever,  must  be  greater  than  the  maximum  gain  per 
stage,  which  is  approximately  20  dB  for  the  three-stage  gyro- 
TWT  amplifier  in  Section  III.  Since  it  is  fairly  straightforward 
to  build  broadband  couplers  and  attenuators  with  a  return  loss 
of  over  10  dB,  this  constraint  is  not  too  demanding. 

It  is  also  assumed  in  developing  the  nonlinear  theory  that 
space  charge  effects  can  be  neglected.  This  assumption  will 
be  valid  for  a  single-stage  device  employing  a  low  value  of 
beam  current.  For  the  harmonic  gyro-TWT  amplifiers  studied 
here,  the  validity  of  this  assumption  remains  to  be  verified 
experimentally,  as  both  high  levels  of  beam  current  and 
multiple  severs  are  required  for  stable  high-power  generation. 

The  validity  of  the  theoretical  prediaion  [17J,  [18]  that 
harmonic  gyro-TWT  amplifiers  can  generate  extremely  high 
power  in  the  millimeter-wave  region  has  been  demonstrated 
through  a  numerical  design  study.  By  limiting  the  length  of 
each  interaction  section,  a  multistage  device  can  be  kept  stable 
and  still  yield  ample  gain.  It  has  also  been  shown  that  a  MIG 
beam  can  be  utilized  with  satisfactory  performance.  A  533 
kW,  35  GHz,  three-stage  second-harmonic  TE;2i  gyro-TWT 
amplifier  with  21.3%  efficiency  and  7.4%  bandwidth  has  been 
simulated  numerically  and  will  be  constructed  at  UCLA.  Fur¬ 
ther  improvement  of  the  performance  of  the  designed  harmonic 
gyro-TWT  amplifier  may  be  expected  if  enhancement  schemes 
such  as  tapering  the  axial  magnetic  field  (28],  [29]  and/or 
waveguide  wall  [5]  are  employed. 
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Abstract — Bragg  reflector  composed  of  periodic  cylindrically 
symmetric  corrugations  can  provide  a  reflection  of  nearly  unity 
through  the  principles  of  constructive  interference,  allowing  the 
formation  of  a  frequency  selective  resonator.  Mode  conversion 
will  occur,  but  can  be  reduced  by  tapering  the  amplitude  of  the 
corrugations.  Reflection  measurements  are  compart  with  theo¬ 
retical  predictions  for  untapered  sinusoidal  and  rectangular  cor¬ 
rugation  as  well  as  a  reflector  tapered  according  to  the  Hamming- 
window  prescription.  Measurements  of  Bragg  resonators  are  also 
presented. 


I.  Introduction 

A  Bragg  resonator  is  the  most  suitable  cavity  structure  for 
providing  frequency  selective  feedback  for  overmoded, 
high-power  oscillators  such  as  the  cyclotron  autoresonance 
maser  (CARM)  [1]  and  the  free  electron  laser  (PEL)  [2].  It 
can  be  obtained  from  corrugated  waveguide  as  shown  in  Fig. 
1.  A  Bragg  resonator  also  allows  good  beam  transport.  Since 
the  corrugation  amplitude  can  be  quite  small,  the  path  for 
the  electron  beam  is  unobstructed.  The  corrugation  scatters 
a  forward  wave  coherently  into  a  backward  wave  if  the  period 
satisfies  the  Bragg  condition 

2k  =  kt,  (1) 
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Fig.  1.  Schematic  of  Bragg  resonator  with  sinusoidal  corrugation. 


mode  number  cannot  change  during  conversion;  therefore  the 
TEii  mode  will  not  couple  to  the  TE21,  TEoi,  or  TMoi 
modes,  but  it  does  to  the  TMn  mode.  If  these  two  modes 
are  close  in  frequency,  then  strong  parasitic  intefcoupling 
can  occur  between  them.  This  places  an  upper  limit  on  the 
corrugation  amplitude  for  a  given  corrugation  geometry. 

The  organization  of  this  paper  is  as  follow.  In  Section  II, 
the  theory  and  measurement  of  Bragg  reflectors  are  given. 
It  is  shown  that  tapering  the  corrugation  can  dramatically 
improve  the  mode  quality.  The  theory  and  measurement  of 
Bragg  resonators  are  presented  in  Section  111,  and  Section  IV 
contains  the  final  discussion  and  conclusions. 


where  k  is  the  axial  propagation  constant  of  the  wave,  kb  = 
2-k/I,  and  I  is  the  corrugation  period.  A  resonator  comprising 
two  Bragg  reflectors  separated  by  a  smooth  tube  as  in  Fig. 
1  can  have  a  high  quality  factor  {Q  >  1000)  for  a  frequency 
satisfying  the  Bragg  condition. 

A  characteristic  of  a  Bragg  reflector  is  that  it  can  also  con¬ 
vert  the  mode  of  a  wave,  especially  if  two  modes  are  close  in 
frequency  for  the  axial  wave  vector  which  satisfies  the  Bragg 
condition.  Our  mode  of  interest  is  the  TEu  mode,  the  lowest 
order  mode  of  cylindrical  waveguide.  The  four  nearest  modes 
are  the  TE21,  TEoi,  TMqi,  and  TMn  modes.  Because  our 
corrugation  geometry  is  cylindrically  symmetric,  the  azimuthal 
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II.  Reflector 

A.  Theory 

The  general  case  of  a  waveguide  whose  cross  section 
varies  smoothly  along  the  axis  has  been  treated  by  Solvmar 
[3],  who  formulated  the  set  of  coupled  equations  for  the 
wave  amplitudes  of  the  forward  and  backward  propagating 
waves,  namely  the  generalized  telegraphist’s  equation.  From 
this  equation,  under  the  assumption  of  a  weak  sinusoidal 
corrugation  one  can  derive  either  a  set  of  coupled  differential 
equations,  where  all  modes  are  assumed  to  be  isolated  from  the 
other  modes,  or  a  more  exact  set  of  multicoupled  differential 
equations,  where  the  modes  are  allowed  to  intercouple.  By 
solving  the  uncoupled  and  coupled  differential  equations,  the 
behavior  of  sinusoidal  corrugation  Bragg  reflectors  can  be 
examined  as  well  as  the  basic  characteristics  of  a  Bragg 
reflector  with  general  corrugation  geometry. 

I )  Single  Mode:  Throughout,  we  will  consider  TEmn  and 
TMmn  modes  whose  fields  satisfy 

/{z,<)  =  (2) 

We  will  assume  a  slow  variation  of  the  amplitude  such  that 
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We  will  also  assume  cylindrical  geometry  with  a  weak  ax- 
isymmetric  sinusoidal  ripple  in  the  waveguide  wall,  so  that 
the  radial  profile  follows 


k{z)  =  a  +  licos{kbz), 

(4) 

where  a  is  the  mean  radius  and  the  corrugation  amplitude 
h  satisfies  h  a.  For  the  single-mode  case,  the  coupled 

equations  for  the  forward  (A"*")  and  backward  (A") 
amplitudes  are 

wave 

dA'*' 

— —  =  — (q  -f-  iA)A'*'  -  iGA  ' 
dz 

(5) 

and 

dA~ 

~ —  =  (q  -1-  iA)A“  +  tGA'*', 

dz 

(6) 

where  A  =  k  —  kb/2  is  the  Bragg  mismatch,  k  =  {(uj/c)^  — 
{^tnn  is  the  axial  wave  vector,  Xmn  represent  the 

zeros  of  the  Bessel  functions  (TM)  and  their  first  derivatives 
(TE),  a  is  the  wall  radius,  a  is  the  ohmic  dissipation  loss 
constant  (it  is  usually  a  good  approximation  to  let  a  =  0), 
and  G  is  the  wave  coupling  coefficient.  The  reflectivity  for  a 
single  seaion  reflector  can  be  obtained  analytically  [4]  for  the 
boundary  conditions 


A-''(0)  =  1 

(7) 

and 

A-(L)  =  0, 

(8) 

where  the  origin  is  taken  to  be  at  one  end  of  the  reflector  of 

length  L.  The  reflectivity  is  found  to  be 


G  sinhyL 

(a  + 1  A)  sinhyL  +  ycoshyL 

(9) 

where  y^  =  +  (q  +  lA)^. 

Fig.  2  shows  the  predicted  reflection  coefficient  from  (9) 
for  a  single-section  reflector  as  a  function  of  Bragg  mismatch 
normalized  to  the  coupling  coefficient.  As  one  can  see  from  the 
figure,  the  bandwidth  is  directly  proportional  to  the  coupling 
coefficient  and  therefore  is  inversely  proportional  to  the  length 
of  the  corrugation  section  for  a  fixed  reflectivity  (fixed  GL). 
Thus,  for  a  narrower  bandwidth  reflector  with  a  specific  desired 
reflectivity,  one  should  reduce  the  corrugation  amplitude  and 
increase  the  length  of  the  reflector  (i.e.,  fix  GL  and  decrease 
G/L  ). 

2)  Multimode:  In  a  given  Bragg  reflector,  all  modes  satisfy 
the  Bragg  condition  at  some  frequency.  In  addition,  the  cor¬ 
rugation  can  couple  a  mode  to  other  modes  if  the  two  waves 
satisfy  the  Bragg  condition 

ki  +  k2  =  kb  (10) 


fZ=|A-(0)/AW  = 


and  the  symmetry  allows  it.  The  most  general  set  of  equations 
for  the  multicoupled  case  is 

^  =  -(a,  +  tA,)A;  -  t  ^  A-  (1 1) 


Fig.  2.  Reflection  coefficient  from  uncoupled  mode  Ihciuy  for  a  section 
of  corrugated  waveguide  with  GL  =  2  versus  Bragg  resonance  mismatch. 


and 

^  =  (a,+iA,)A-+t£Gj.A^’  (12) 

• 

where  the  subscript  refers  to  the  various  transverse  modes  and 
Gij  is  the  wave  coupling  coefficient  tensor. 

Since  our  mode  of  interest  is  the  TEn  mode,  we  will 
confine  the  rest  of  the  analysis  to  just  the  TEn  and  TMn 
modes.  Then  the  general  equations  for  the  multicoupled  case 
reduce  to  four  first-order  linear  equations  (H  representing 
TEii  and  E  representing  TMn): 


dA 


dz 


:= -(off  +  iAH)A^  -  iGffAJji  -  iGxAjr  (13) 


dA 


dz 

dA+ 


H  _ 


=  {ofi  -I-  zA//)Af^  +  iGiiA'fj  +  iGxA^ 


and 


dz 

dA 


—  =  —{oE  -f  lAfrjA^  —  iGsA^  -  iGxAjj 


(14) 

(15) 


dz 


B  _ 


=  {ocE  +  iAE)A^  +  iGeA^  +  iGxA'^i  (16) 


where  again  Ah  and  A^  represent  the  Bragg  mismatches 
for  the  TEii  and  TMn  modes,  respectively,  and  Gu,  Ge, 
and  Gx  am  the  wave  coupling  coefficients  for  TEn  “  TEn, 
TMn  -  TMn,  and  TEn  -  TMn  conversion,  respectively. 
The  general  coupling  coefficients  for  sinusoidal  corrugation 
are  given  by  [5] 

h  -1-  mmaHu>yc‘^  -F  kk‘)) 

"  2a3  (fc):‘)t/2((x^„  -  m2)(x^„.  - 

for  the  coupling  of  TEmn  with  TEm  n'  modes, 

h  (w^/c^  +  kk‘) 

2^  {kk'yi^ 

for  the  coupling  of  TMmn  with  TM^  n-  modes,  and 


Ge  = 


(18) 


_  /l  m(uj/c)  (•'  +  k  ) 


(19) 


for  the  coupling  of  TEmn  with  TMm  n  modes.  If  the  oper¬ 
ating  frequency  is  much  greater  than  the  cutoff  frequency,  the 
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coupling  coefticient.s  can  be  siniplilied  as 


TABLI  I 

COKKl'I.AnON  Kinjci'oi'  pAHAMm_RS 


t;,.  -  II  12  -  -  (20) 

(i  (■ 

i'l.n  mode. 

h 

(21) 

C'jso 

:so  GHz 

UK.  C.He 
(bl 

to  b  GHz 
(c) 

Period 

0.f>03  mm 

1.575  cm 

1 .575  cm 

Mean  radius 

3.525  mm 

1  cm 

1.809  cm 

il  ( 

r.\u;  mode,  and 

Corrugaiion 

amplKudc 

0  02b  mm 

0  cm 

0.203  cm 

h 

(22) 

Corrugation  length 

181)^  cm 

10 

40. 8o  cm 

1  2 

31.50  cm 

? 

a  r 

Corrugation  profile 

constant 

constant 

constani 

intercoupling  resonance  From  (17)-(19)  or  (20)-(22), 

Corrugation  gcomcirv 

sinusoidal 

sinusoidal 

reaanguiar 

to.'  the 


one  can  observe  that  the  coupling  strength  is  proportional  to 
the  normalized  corrugation  amplitude  li/a. 

Athough  we  cannot  obtain  a  simple  closed-form  solu¬ 
tion  of  the  four  coupled  equations,  a  numerical  solution  can 
be  obtained  for  the  wave  amplitudes  by  using  the  eigen- 
value/eigenvector  method.  Tire  four  coupled  equations  can  be 
expressed  in  vector  form  as  shown  in  (23),  at  the  bottom  of 
the  page.  Letting  the  amplitudes  satisfy  A  <x  where  A 
is  a  complex  constant,  one  can  express  the  coupled  vector 
differential  equations  as  a  compact  4x4  matrix  eigenvalue 
problem  given  by 


for  TEii(TMii)  input. 

Using  this  analytical  theory,  the  reffectivity  characteristics 
of  sinusoidal  corrugation  have  been  examined.  The  reflectivity 
for  TEii  and  TMji  modes  is  given  by 

Rh  -  |>lw(0)i'  (30) 


and 


(31) 


(A4-AJ)e  =  0,  (24) 

where  M  is  the  4x4  matrix  on  the  RHS  of  (23),  A  is  an 
eigenvalue,  e  is  the  corresponding  eigenvector,  J  is  the  4x4 
identity  matrix,  and  0  is  the  zero  vector.  One  can  recognize  that 
(24)  is  in  the  familiar  eigenform  and  therefore  the  eigenvalues 
and  eigenvectors  can  straightforwardly  be  obtained. 

Once  the  eigenvalues  and  eigenvectors  have  been  found, 
the  amplitude  vector  A  at  any  position  of  the  corrugation 
section  can  be  expressed  as  a  linear  combination  of  the  four 
eigenvectors  as 

A  =  (25) 

1=1 

where  K,  arc  constants  to  be  found  by  satisfying  the  boundary 
condition.  The  equations  have  been  solved  for  the  boundary 
conditions 


A*(0]  =  1(0) 

(26) 

AlJ(O)  =  0(1) 

(27) 

AJf(L)  =  0(0) 

(28) 

Al{L)  =  0(0) 

(29) 

Similarly,  the  transmission  through  the  reflector  is  given  by 

'Rh  —  Mh(^)|  (3^) 

and 

Tf  =  |A^(L)f.  (33) 

Fig,  3(a)  shows  the  reflection  characteristics  for  TEn  input 
from  a  sinusoidal  corrugation  Bragg  reflector  with  the  param¬ 
eters  listed  in  Table  I  (case  (a)).  This  reflector  was  built  at  the 
Lawrence  Livermore  National  Laboratory  (LLNL)  for  a  250 
GHz  CARM  oscillator  experiment  [6],  (7].  If  the  amplitude 
is  increased  from  the  value  in  Fig.  3(a),  the  reflectivity  in  the 
TEii  mode  actually  decreases  because  the  resonances  begin 
to  overlap  and  interfere  with  one  another.  Notice  in  Fig.  3(b) 
that  intercoupling  destroys  the  performance  of  the  reflector  for 
a  large  corrugation  amplitude.  The  reflectivity  for  TMjj  mode 
input  for  the  same  reflector  is  shown  in  Fig.  3(c).  Notice  that 
the  TEii  -  TMii  intercoupling  resonanc'  occurs  at  the  same 
frequency  in  Fig.  3(c)  as  in  Fig.  3(a),  as  would  be  expiected 
from  reciprocity.  It  is  observed  in  Fig.  3(d)  that  the  mode 
purity  of  the  TMn  mode  is  also  degraded  as  a  result  of 
overlap  with  the  TEn  -  TMu  intercoupling  resonance  for 
a  large  corrugation  amplitude. 

The  maximum  allowable  corrugation  amplitude  for  a 
high  mode  purity  reflector  is  determined  by  the  need 


JZ 


'^H  \  1 

(  -{qh  +  »2)sw) 

0 

-tGx 

1  - 

iGn 

{a If  -f  jAh) 

iG.x 

0 

1 

0 

-rG.v 

-(«■£  +  *Ae) 

—iG  £• 
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V  tGx 
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Fig.  3.  Dependence  on  frequency  of  reflection  in  TEu  (solid  line)  and  TMn  (dashed  line)  modes  from  sinusoidal  corrugation 
reflector  with  parameters  otherwise  corresponding  to  those  in  Table  I,  case  (a),  with  corrugation  amplitudes  of  (a)  0.026  mm  and 
(b)  0.039  mm  for  TEn  mode  input  and  (c)  0.026  mm  and  (d)  0.039  mm  for  TMn  mode  input 


to  minimize  the  overlap  of  the  nearest  intercoupling 
resonance  with  the  desired  Bragg  resonance.  The  fre¬ 
quency  separation  between  the  two  modes  must  be  at 
least  greater  than  the  sum  of  the  half-widths  of  each 
resonance.  For  the  case  of  interest  and  using  Afc  = 
3G  as  the  full  width  (see  Fig.  2),  small  mode  overlap 
requires 

-^H>  (3/2)c(|Gh|  +  |Gx|),  (34) 

where  wfj  and  ujx  are  the  center  frequencies  of  the  pure 
TEii  and  TEn  -TMu  resonances,  respectively.  Multiplying 
both  sides  by  (wx  +  wh)  —  2ujf/,  we  obtain  the  condi¬ 
tion 

ori- -a;?,  >3a;„c(|G//|  +  |Gx|).  (35) 


By  using  the  dependence  of  frequency  on  the  cutoff 
frequency  and  propagation  constant,  it  can  be  shown 
that 


4 


where  cjce  and  ui^h  are  the  cutoff  frequency  of  the  TMn  and 
TE,i  modes,  respectively.  Substituting  (20),  (22)  and  (36) 


into  (35)  and  expressing  u-ce  in  terms  of  u-ch,  we  find  that 
the  maximum  allowable  corrugation  amplitude  is  determined 
by  the  condition 

h/a  <  0.78(uJc///u;)^  (37) 

It  is  evident  that  the  further  the  device  is  operated  above 
cutoff,  the  smaller  the  maximum  allowable  corrugation  am¬ 
plitude. 

Bragg  reflectors  are  designed  by  setting  the  corru¬ 
gation  amplitude  to  a  value  at  which  the  resonances 
are  well  separated  and  then  choosing  the  length  of  the 
reflector  to  obtain  the  desired  reflection  Fig.  4  shows 
the  depiendence  on  length  of  Bragg  reflection  for  the 
250  GHz  LLNL  Bragg  reflector  as  predicted  by  the 
coupled  mode  theory.  Reducing  the  corrugation  amplitude 
further  would  yield  a  smaller  fraction  in  the  TMu 
mode. 

3)  Scattering  Matrix  Code:  Since  the  frequency  selectivity 
of  a  Bragg  reflector  depends  only  on  the  axial  periodicity  of 
the  waveguide  wall  perturbation,  the  corrugation  geometry, 
in  general,  can  be  chosen  arbitrarily.  However,  the  previous 
analytical  theory,  which  assumed  a  slow  and  weakly  varying 
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Fig.  4.  Dependence  on  corrugation  length  of  reflection  at  (a)  249  GHi  and 
(b)  250  GHz  in  TEn  (solid  lines)  and  TMn  modes  (dashed  lines)  from 
sinusoidal  corrugation  reflector  (Table  I,  case  (a)).  The  injected  wave  is  in 
TEii  mode. 

corrugation,  is  not  truly  appropriate  for  a  Bragg  reflector  with 
step  discontinuities.  To  study  the  reflectivity  characteristics  of 
rectangular  cornjgation  reflectors,  in  lieu  of  determining  the 
correct  coupling  coefficients,  a  numerical  scattering  matrix 
simulation  is  used.  The  scattering  matrix  method  (8],  (9]  is 
a  multimode  formulation  where  many  propagating  modes  are 
included  as  well  as  the  evanescent  modes.  Inclusion  of  the 
evanescent  modes  is  especially  important  when  they  propagate 
as  real  waves  in  the  corrugation  vanes. 

In  this  approach,  the  waveguide  is  decomposed  into  a  dis- 
Crete  number  of  modules,  corresponding  to  the  number  (minus 
one)  of  radial  discontinuities.  Each  consists  of  an  abrupt  junc¬ 
tion  of  two  circular  waveguides  with  different  diameters.  If  the 
reflector  is  externally  excited  in  the  fundamental  TEu  mode, 
the  complete  set  of  eigenfunctions  for  each  module  is  repre¬ 
sented  by  the  forward  and  backward  traveling  and  evanescent 
cylindrical  waveguide  TEi„  and  TMi„  modes.  Each  TEjn 
or  TMj,.  wave  incident  on  a  discontinuity  excites  an  infinite 
sum  of  transmitted  and  reflected  TEjn  and  TMi„  waves 
as  determined  by  the  coupling  coefficients  at  the  interface. 
Truncation  of  this  system  for  the  numerical  evaluation  leads 
to  the  individual  scattering  matrix  of  each  corrugation  module. 

The  complete  scattering  matrix  for  the  entire  reflector 
is  determined  by  repetitively  cascading  the  modules.  Finite 
conductivity  is  included  in  the  computer  code  to  model 
the  performance  of  the  Bragg  reflector  accurately.  Since 
the  fields  in  each  module  are  expressed  by  a  sum  of  the 
eigenfunctions  of  smooth  waveguide  with  a  constant  radius, 
rectangular  corrugation  is  the  most  convenient  geometry  for 
this  calculation.  However,  this  method  can  be  extended  to 
the  simulation  of  any  periodic  form  by  dividing  a  corrugation 


Frequency  (GHz) 

Fig.  5.  Frequency  dependence  of  insertion  loss  through  sinusoidal  corruga¬ 
tion  reflector  (Table  1.  case  (a))  of  TEn  mode  from  two-mode  theory  (solid 
line)  and  measurement  (circles) 

with  an  arbitrary  shape  into  many  small  steps  of  rectangular 
discontinuities  to  approximate  the  particular  geometry. 

B.  Measurement 

«  , 

Three  sets  of  constant-corrugation-amplitude  Bragg  reflec¬ 
tors  have  been  fabricated  and  tested.  Two  sets  of  Bragg 
reflectors  designed  for  10.6  GHz  with  ujjuJc  =  2.2  were  built 
and  tested  at  UCLA  and  another  set  designed  for  250  GHz  with 
ui/wc  =  10  was  built  and  tested  at  LLNL.  All  three  sets  were 
electroformed  with  copper.  The  inner  surfaces  of  the  UCLA 
reflectors  were  plated  with  gold  to  avoid  oxidation  during  the 
test  measurement  phase. 

The  insertion  loss  of  the  250  GHz  sinusoidal  corrugation 
reflector  specified  by  Table  1  (case  (a))  was  measured  by 
attaching  a  tunable  source  to  the  reflector  through  a  down- 
taper  to  fundamental  rectangular  waveguide.  The  ether  end 
of  the  corrugation  section  was  left  open.  The  transmitted 
signal  was  picked  up  by  a  horn  and  measured  with  a  large 
dynamic  range  heterodyne  receiver.  The  measured  insertion 
loss  is  shown  in  Fig.  5  and  should  be  compared  with  the 
predicted  transmissivity  from  the  analytical  coupled  mode 
theory  also  shown  in  Fig.  5.  The  resonance  at  248.8  CHz 
is  due  to  pure  TEn  coupling  and  the  resonance  at  250.7  GHz 
is  from  TEn  -  TMn  intercoupling.  The  measurements  are  in 
good  agreement  with  the  theoretical  predictions  including  the 
small  gap  between  the  desired  TEn  —  TEn  resonance  and 
the  destructive  TEn  -  TMn  intercoupling  resonance.  Even 
though  the  reflector  was  designed  for  a  large  overmoded  value 
of  w/w„  =  10,  the  mode  overlap  is  at  a  tolerable  level  because 
of  the  corrugation’s  small  amplitude. 

The  return  loss  of  each  UCLA  Bragg  reflector  was  measured 
with  an  HP  8720A  automated  vector  network  analyzer.  The 
measurement  setup  is  shown  is  Fig.  6.  A  down-taper  is 
attached  to  one  end,  followed  by  a  circular  to  TEm  rectangular 
waveguide  transformer,  while  the  other  end  is  left  open.  For 
frequencies  within  the  design  bandwidth,  the  converted  TMn 
mode  will  be  cut  off  in  the  down-taper.  Neglecting  interference 
effects,  the  total  rcflecied  power  picked  up  by  the  receiver  can 
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Fig  (i  Schematic  diagram  of  the  UCLA  measurement  setup  to  determine 
the  Bragg  reflection  coefficient  (Si,) 

be  expressed  as 

OC 

■^TOTAL  =  /?//  +  Rx  (38) 

1=0 

where  Rfj  and  R^  are  the  single  pass  reflection  from  TEn 
and  TMii  modes,  respectively,  Rx  is  the  mode  conversion 
between  TEn  and  TMn  modes,  and  firoTAL  is  the  the  total 
reflected  power  into  the  TEn  mode. 

7(a)  shows  the  predicted  reflectivity  from  the  analytical 
coupled  mode  theory  for  the  sinusoidal  corrugation  reflector 
specified  in  Table  I  (case  (b)),  TTie  measurement  of  the  return 
loss  (5n)  for  this  reflector  is  shown  in  Fig.  7(b)  and  should 
be  compared  with  the  solid  line  in  Fig.  7(a),  which  includes 
the  effect  of  multiple  reflection  (eq.  (38)).  The  agreement  of 
the  measurement  with  the  theoretical  prediction  is  fairly  good 
except  that  the  experimental  bandwidth  is  25%  broader  than 
the  theoretical  prediction. 

The  Bragg  reflector’s  corrugation  need  not  be  sinusoidal.  A 
reflector  was  also  built  and  tested  with  a  rectangular  profile. 
Fig.  8(a)  shows  the  theoretical  reflectivity  of  the  rectangular 
corrugation  reflector  specified  by  Table  I  (case  (c))  as  predicted 
by  the  scattering  matrix  simulation  code.  Once  again,  the 
solid  curve  includes  the  effect  of  multiple  reflection.  The 
measurements  of  this  reflector  are  shown  in  Fig.  8(b).  Notice 
that  the  10.5-12.5  GHz  TEh  — TMh  intercoupling  resonance 
interferes  with  the  desired  TEjj  —  TEn  coupling  resonance  at 
9.5-1 1  GHz,  as  predicted  by  the  scattering  matrix  simulation. 
The  experimental  results  are  in  very  good  agreement  with  the 
scattering  matrix  simulation  results,  including  the  bandwidth  of 
the  two  resonances.  However,  the  analytical  theory  (Fig.  8(c)) 
using  the  larger  (x4/7r)  equivalent  fundamental  sinusoidal 
corrugation  amplitude  obtained  from  Fourier  decomposition, 
which  does  not  include  the  cutoff  modes  in  the  analysis,  failed 
to  predict  the  extended  bandwidth  below  the  cutoff  frequency 
for  the  TMji  mode  (10.1  GHz)  as  well  as  the  stronger 
coupling  characteristics  of  the  rectangular  corrugation. 

C.  Hamming-Window 

From  the  simulation  results  for  the  previous  constant- 
amplitude  corrugation  reflectors,  wc  have  observed  that  the 
side  lobes  of  the  intercoupling  resonance  destroys  the  mode 
purity  of  the  pure  TEn  -  TEn  resonance.  The  problem 
of  undesirable  side  lobes  can  be  reduced  by  tapering  the 
corrugation  amplitude  profile  as  pre.scribcd  by  a  Hamming- 
window  distribution  (10],  a  common  technique  in  filter  theory 
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Fig.  7.  Frequency  dependence  of  refleciion  from  sinusoidal  corrugation 
reflector  (Table  1,  case  (b))  (a)  in  TE,,  (thin  line)  and  TMn  (dashed  line) 
modes  from  two-mode  theory  including  total  theoretical  TEn  reflection  (dark 
line)  from  multiple  reflections  and  (b)  from  measurement.  The  input  wave  is 
in  the  TEn  mode 

to  suppress  the  side  lobes  of  a  filter’s  response.  The  envelope 
of  the  Hamming-window  amplitude  distribution  is  a  modified 
raised  cosine  distribution  given  by 

h{z)  =  /io(0.54  -0.46cos(27r2/L)).  (39) 

where  h{z)  is  the  amplitude  of  the  corrugation  at  position  c 
(distance  being  measured  from  the  beginning  of  the  reflector 
of  length  L)  and  /lo  is  the  peak  amplitude  of  the  corrugation. 

We  have  tested  the  reflectivity  of  two  Hamming-window 
Bragg  reflectors  with  rectangular  corrugation  for  the  gun  and 
collector  ends  of  the  cavity  for  a  planned  400  kV  CARM 
oscillator  experiment  at  10.6  GHz  R-j/wc  =  2.2)  [llj.  They 
were  also  electroformed  with  copper.  We  have  chosen  rectan¬ 
gular  corrugation  for  our  Hamming-window  reflector  because 
the  rectangular  geometry  is  mechanically  easier  to  fabricate 
than  the  curved  surfaces  of  sinusoidal  corrugation  and  the 
change  from  untapered  to  tapered  amplitude  profile  in  the 
scattering  matrix  formulation  is  a  relatively  straightforward 
modification.  The  corrugation  orientation  of  our  Hamming- 
window  reflectors  is  positive  (Il(z)  >  0)  to  avoid  creating  a 
high-(>  cavity  at  the  cutoff  frequency.  The  positive  corrugation 
also  provides  a  larger  open  passage  for  the  electron  beam  to 
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Fig  9.  Frequency  dependence  of  reflection  from  Hamming-wmdow  gun-side 
reflector  (Table  II)  (a)  In  TEn  (solid  line)  and  TM|i  (dashed  line)  modes 
from  scalteriiig  matrix  simulation  and  (b)  from  measurement.  The  input  wave 
is  in  the  TEn  mode 


(c) 

Fig  8  Frequency  dependence  of  reflection  from  rectangular  corrugation 
reflector  (Table  1.  case  (c))  (a)  in  TEn  (thin  line)  and  TMn  (dashed  line) 
modes  from  scattering  matrix  simulation  including  total  TE|  |  reflection  (dar)- 
line)  from  multiple  reflections,  (b)  from  measurement,  and  (c)  from  two-mode 
theory  The  input  wave  is  in  TEn  mode. 


0.10%  in  the  TMn  mode,  while  the  output  reflector  transmits 
3.41%  in  the  TEn  mode  and  0.11%  in  the  TMn  mode. 

III.  Resonator 


pass  through.  The  positive  corrugation  amplitude  is  defined 
with  respect  to  the  minimum  inner  wall  radius. 

Figs.  9(a)  and  10(a)  show  the  scattering  matrix  simulation 
results  for  the  gun-side  and  output  Bragg  reflectors,  respec¬ 
tively,  described  by  Table  II.  Ohmic  loss  is  included.  The 
simulation  found  a  peak  reflectivity  of  99.4%  for  the  gun-side 
reflector  and  96.0%  for  the  output  reflector.  For  an  injected 
TEji  wave,  the  simulation  also  found  that  the  reflected  wave 
at  10.6  GFlz  from  both  the  gun-side  and  the  output  reflector  is 
99.9%  TEii-  The  measured  5iivalues  for  the  corresponding 
reflectors  are  shown  in  Figs.  9(b)  and  10(b)-  As  shown  in  the 
figures,  the  results  for  the  reflectivity  and  bandwidth  are  seen 
to  be  in  excellent  agreement  with  theory.  As  predicted  by  the 
simulation,  the  frequency  dependence  of  the  reflection  for  both 
tapered  reflectors  is  very  smooth  compared  with  the  untapered 
sinusoidal  corrugation  measurement  (Fig.  7(b)). 

The  transmission  characteristics  of  the  Hamming-window 
reflectors  are  shown  in  Fig.  11.  The  results  from  the  scattering 
matrix  simulation  are  shown  in  parts  (a)  and  (b)  of  the  figure 
for  the  gun-side  and  output  reflectors,  respectively.  At  the 
design  frequency  of  10.6  GHz,  the  gun-side  reflector  transmits 
only  0.08.S%  of  the  input  TEn  wave  in  the  Ti:,,  mode  and 


A.  Theory 

A  Bragg  resonator  cavity  can  be  formed  by  separating  two 
corrugation  sections  by  a  smooth  section  of  length  Lr.  The 
diffractive  quality  factor  for  the  cavity  can  be  approximated 
by  the  well-known  expression  (1) 

Q^[u>lr.?LRl[k[\- R^R2)\■ 

where  R\  and  /?2  are  the  reflectivities  (in  amplitude)  of  the  two 
corrugated  sections.  However,  the  above  result  is  only  valid  for 

GLr-»\ 

~  1  (42) 

/?2~1.  (43) 

A  better  prediction  for  Q,  including  intercoupling  as  well  as 
the  energy  resonating  inside  the  reflectors  themselves,  can  be 
obtained  by  using  the  previous  eigenvalue/eigenvecior  method 
(eq.  (23))  for  the  sinusoidal  corrugation  reflectors.  In  this 
analysis,  the  wave  is  injected  through  one  of  the  two  reflectors 
In  the  resonating  section,  the  uncoupled  amplitudes  satisfy 

dA'jd:  =  Tin,  ■+-  (44) 

or  (4.S),  shown  at  the  hottoin  I'l  the  next  page. 
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Frequency  (GHz) 


(b) 

Fig.  10.  Frequency  dependence  of  reflection  from  Hamming-window  output 
reflector  (Table  II)  (a)  in  TEn  (solid  line)  and  TMn  (dashed  line)  modes 
from  scattering  matrix  simulation  and  (b)  from  measurement.  The  input  wave 
is  in  the  TE|  i  mode. 


TABLE  II 

Parameters  of  Hammino  Window  Brago  Reflectors 
Fabricated  for  10.6  GHz  CARM  Osollator  Experiment 


Gun-Side 

Output 

Peri(xl 

1.500  cm 

IJOCcm 

Mean  radius 

1 .809  cm 

1.809  cm 

Coirugation 

amplitude 

0.400  cm 

0.325  cm 

Conugalion  length 

36.75  cm 

36.75  cm 

2.2 

2.2 

Corrugation 

orientation 

positive 

positive 

Corrugation  profile 

tapered 

tapered 

Corrugation  geometry 

rectangular 

rectangular 

The  matrix  M  for  the  smooth  section  is  therefore  composed 
of  only  the  propagating  constants  for  the  TEu  and  TMn 


(b) 

Fig.  11.  Frequency  dependence  of  TEu  (solid  lines)  and  TMn  (dotted 
lines)  mode  transmission  through  Hamming-window  reflectors  (Table  II)  from 
scattering  matrix  simulation  for  (a)  gun-side  reflector  and  (b)  output  reflector. 
The  input  wave  is  in  the  TEn  mode 

modes.  There  are  four  eigenvalues  and  eigenvectors  in  the 
resonator  section  and  each  of  the  reflectors  and  therefore 
a  total  of  12  unknown  coefficients.  A  unique  solution  for 
the  wave  amplitudes  can  be  obtained  when  12  boundary 
conditions  are  specified.  In  addition  to  the  four  input  and 
output  boundary  conditions,  the  eight  additional  boundary 
conditions  are  realized  by  the  continuity  of  the  waves  at  the 
two  ends  of  the  resonating  section. 

Fig.  12  shows  the  predicted  frequency  response  of  energy 
resonating  in  a  Bragg  cavity  with  two  Bragg  reflectors  de¬ 
scribed  by  Fig.  3(a)  and  separated  by  a  smooth  section  of 
length  Lfi  —  b  cm.  The  calculation  includes  an  ohmic  loss 
of  twice  the  theoretical  value  in  both  the  reflectors  and  the 
central  resonator  section.  The  modes  at  249.0  GHz  and  250.2 
GHz  are  predominantly  TEu,  while  the  mode  at  251.0  GHz 
is  approximately  a  50%  mixture  of  TEu  and  TMu-  The  Q 
of  a  mode  in  the  resonator  can  be  calculated  as 

<3=1^.  (46) 

where  /o  is  a  resonant  frequency  and  A/  is  the  half-power 
width.  This  gives  Q  values  of  1 1  000  and  4000  for  the  mcxles 
at  249  GHz  and  250.2  GHz,  respectively. 


dz 


Ah  ' 

/  -{an  -f  ikn) 

0 

0 

0 

-  1  ° 

(oh  +  rkf]) 

0 

0 

0 

0 

-{<»£:  +  lAf) 

0 

\aiJ 

\  0 

0 

0  (, 

(45) 


..  UKAOC)  Kh^H;CTORS 


401 


100 

75 

50 

O 

o 

25 

a- 

0 

248 


249  250  251  252 

Frequency  (GHz) 


Frequency  (GHz) 


Fig.  12.  Frequency  response  of  resonator  with  two  sinusoidal  corrugation 
reflectors  described  in  Table  I  (case  (a))  for  Lr  =  5.0  cm  (ohmte  loss 
included). 


Fig  13.  Spattal  variation  of  the  mode  from  Ftg.  12  at  249.0  GHz  for  Z.r  * 
5  cm.  The  dashed  line  shows  the  envelope  of  the  energy  profile  and  the  solid 
line  shows  the  energy  profile  itself. 


The  Spatial  profile  of  the  central  mode  (/=  249.0  GHz)  from 
the  resonator  described  in  Fig.  12  is  shown  in  Fig.  13.  It  is 
evident  that  a  significant  amount  of  energy  leaks  into  each 
reflector  from  the  cavity  region. 


B.  Measurement 

Two  Bragg  resonators,  one  with  the  sinusoidal  corrugation 
reflectors  described  in  Table  I  (case  (b))  and  the  other  with  the 
two  dissimilar  Hamming-window  Bragg  reflectors  described  in 
Table  II,  were  formed  by  separating  two  reflector  sections  by 
a  smooth  uncorrugated  aluminum  tube.  The  inner  radius  of  the 
smooth  central  section  is  equal  to  the  minimum  inner  radius  of 
the  corrugation  for  both  resonators.  The  transmission  method 
has  been  used  to  measure  the  resonant  frequency  and  Q  of 
the  modes  in  both  cavities.  The  signal  was  coupled  in  and  out 
through  piorts  on  the  sidewall  of  the  smooth  resonator  section. 

The  transmission  method  measures  the  loaded  Q  of  the 
resonating  system,  given  by 


Ql  Qsragg  Qohmic  Qext 

In  order  to  directly  measure  the  quality  factor  Q  of  the  Bragg 
cavity,  Q„t  must  be  large  compared  with  (?Br»gg  The  value 
of  Qext  can  be  adjusted  by  decreasing  the  coupling  strength 
of  the  sidewall  coupling  port  such  that 

Qtxt  ^  QBr»gg-  (48) 


Fig  14  Theorelical  dependence  of  a  sinusoidal  Bragg  re.sonaior's 
quality  factor  on  frequency  from  measurement  (open  circles),  the 
eigenvalue/etgenvector  method  (solid  line),  and  eq  (40)  (dashed  line)  The 
reflectors  are  described  in  Table  I  (case  (b))  and  the  variable  smrxith  section's 
length  is  approximately  23  cm  (ohmic  loss  included) 


Then,  the  effective  loaded  Q  of  the  resonator,  neglecting  Qexi< 
is  given  by 


^Bragg  *?Bragg  Qohmic 

Both  the  experimental  and  theorelical  values  of  Qsragg 
been  found  over  the  bandwidth  of  the  sinusoidal  corrugation 
Bragg  reflector  by  changing  the  length  of  the  smooth  section 
of  the  resonator,  thereby  changing  the  resonant  frequency 
of  the  mode.  The  theoretical  values  are  calculated  for  the 
TEi.1,15  axial  mode.  The  measurements  of  QQ,^gf.  are  shown 
in  Fig.  14  along  with  the  theoretical  predictions  from  the 
eigenvalue/eigenvector  method,  which  includes  ohmic  losses 
equal  to  twice  the  theoretical  value,  and  the  qualitative  approx¬ 
imation  (eq.  (40))  using  the  theoretical  values  of  li  from  Fig. 
7(a).  The  experimental  bandwidth  (or  <?BragK  broader  than 
the  theoretical  prediction,  as  was  the  experimental  bandwidth 
for  reflectivity. 

The  measured  Qsiagg  resonator  with  Hamming- 

window  Bragg  reflectors  is  shown  in  Fig.  15  along  with 
the  qualitative  approximation  (eq.  (40))  using  the  reflectivity 
predicted  by  the  scattering  matrix  code.  The  measured  Qeragg 
over  the  bandwidth  is  in  very  good  agreement  with  the 
theoretical  prediction.  Considering  the  sensitivity  of  (40)  on 
1/(1  —  |/?i||/?2l)  ^or  |/?i||/?2|  ~  0.95,  the  small  discrepancy- 
in  the  Q  values  is  quite  reasonable, 

IV.  Conclusions 

Measurements  of  the  reflectivity  characteristics  of  sinusoidal 
corrugation  and  rectangular  corrugation,  both  tapered  and 
untapered,  have  been  compared  favorably  with  theory.  It  was 
shown  that  both  sinusoidal  and  rectangular  corrugation  can 
be  used  as  a  Bragg  reflector.  Tapered  corrugation  was  shown 
to  reduce  mode  conversion  in  the  reflector  and  thus  yield  a 
reflector  with  high  mode  purity. 

Q  measurements  of  Bragg  resonators  with  sinusoidal  cor¬ 
rugation  reflectors  and  Hamming-window  reflectors  were  also 
presented.  It  has  been  shown  that  Bragg  reflectors  can  provide 
a  high-d^,  frequency  selective  resonator  with  good  mode 
purity. 
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Frequency  (GHz) 


Fig  IS.  Theorelical  dependence  of  a  Hamming-window  Bragg  resonaior’s 
quality  factor  on  frequency  from  measurement  (open  circles)  and  eq  (40) 
(dashed  line).  The  two  dissimilar  reflectors  are  described  in  Table  II  and 
the  variable  smooth  section's  length  is  approximately  24  cm  (ohmic  loss 
included). 
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